Methods of specifically labeling nucleic acids using crispr/cas

ABSTRACT

Provided herein are methods of detecting a target nucleic acid sequence. In one embodiment, the method includes contacting genomic DNA with a guide RNA having a portion complementary to the target sequence in the genomic DNA and with Cas9 nickase to produce a single-strand break in the genomic DNA at a specific location adjacent to the target sequence. The method further includes contacting the nicked DNA with a polymerase and fluorescently labeled nucleotide. The fluorescently labeled nucleotide is incorporated into the nicked DNA at the specific location and the target nucleic acid sequence is detected via fluorescent label.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims the benefit under 35 USC 119(e) of U.S. Provisional Patent Application No. 62/410,322, filed Oct. 19, 2016, which is incorporated by reference herein.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

This invention was made with government support under R01 HG005946 and CA177395 awarded by National Institutes of Health. The government has certain rights in the invention.

INCORPORATION-BY-REFERENCE OF MATERIAL SUBMITTED IN ELECTRONIC FORM

Applicant hereby incorporates by reference the Sequence Listing material filed in electronic form herewith. This file is labeled “DRX 15-1773US_ST25.txt”.

BACKGROUND OF THE INVENTION

Two of the major challenges in genome analysis are de novo genome sequence assembly based on ‘short read’ shotgun sequencing and structural variation analysis. Several approaches and combinations of different approaches have been attempted to meet these challenges. The most widely adopted strategy relies on deep sequencing of shotgun libraries and sequencing of mate-pair libraries, which increases the sequence contiguity of short-read sequencing (See, Siegel, A. F., et al. (2000) “Modeling the feasibility of whole genome shotgun sequencing using a pairwise end strategy.” Genomics 68(3): 237-246). The paired sequencing approach includes conventional mate-pair libraries, labor-intensive fosmid or BAC clone libraries (See Gnerre, S., et al. (2011) “High-quality draft assemblies of mammalian genomes from massively parallel sequence data.” Proceedings of the National Academy of Sciences of the United States of America 108(4): 1513-1518), Hi-C read-pairs for chromosome-scale scaffolding (See Burton, J. N., et al. (2013) “Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions.” Nature Biotechnology 31(12): 1119-1125) and transposase-mediated libraries (See Adey, A., et al. (2014) “In vitro, long-range sequence information for de novo genome assembly via transposase contiguity.” Genome Research 24(12): 2041-2049). Another approach relies on the stochastic separation of corresponding genomic or polymerase chain reaction (PCR) fragments into physically distinct pools followed by subsequent fragmentation to generate shorter sequencing templates (See, Kaper, F., et al. (2013). “Whole-genome haplotyping by dilution, amplification, and sequencing.” Proceedings of the National Academy of Sciences of the United States of America 110(14): 5552-5557; Kuleshov, V., et al. (2014) “Whole-genome haplotyping using long reads and statistical methods.” Nature Biotechnology 32(3): 261-266; Peters, B. A., et al. (2012) “Accurate whole-genome sequencing and haplotyping from 10 to 20 human cells.” Nature 487(7406): 190-195; and Voskoboynik, A., et al. (2013) “The genome sequence of the colonial chordate, Botryllus schlosseri.” Elife 2(e00569)). With appropriate high-throughput reaction handling and barcoding, this strategy reduces the complexity, and thus can improve the quality, of assemblies. Longer-read sequencing technologies such as PacBio®'s SMRT and Oxford Nanopore sequencing promise to eventually further improve assembly contiguity. For example, SMRT sequencing has been successfully applied to closing some gaps and detecting some structural variations in the human reference genome (For example, See Chaisson, M. J. P., et al. (2015) “Resolving the complexity of the human genome using single-molecule sequencing.” Nature 517(7536): 608-611). However, their high error rate, low throughput and high cost have thus far prevented widespread adoption.

None of the aforementioned approaches, however, adequately address the problems of long-range de novo assembly contiguity and validation, sequence mis-assembly in complex segmentally duplicated and repetitive regions, and structural variant detection and delineation. Whole genome mapping technologies have been developed for these purposes as complementary tools to provide scaffolds for genome assembly and structural variation analysis. Optical mapping, pioneered by David Schwartz and colleagues has been used to construct restriction maps for various genomes and has proven to be very useful in providing scaffolds for shotgun sequence assembly and detection of structural variations (See, Samad, A., et al. (1995) “Optical Mapping—A novel, single-molecule approach to genomic analysis.” Genome Research 5(1): 1-4; and Teague, B., et al. (2010) “High-resolution human genome structure by single-molecule analysis.” Proceedings of the National Academy of Sciences of the United States of America 107(24): 10848-10853). Furthermore, Ming Xiao and colleagues developed a highly-automated whole genome mapping in a nanochannel array (Hastie, A. R., et al. (2013). “Rapid Genome Mapping in Nanochannel Arrays for Highly Complete and Accurate De Novo Sequence Assembly of the Complex Aegilops tauschii Genome.” Plos One 8(2); Lam, E. T., et al. (2012) “Genome mapping on nanochannel arrays for structural variation analysis and sequence assembly.” Nature Biotechnology 30(8): 771-776 and US 2016/0168621 A1. Each of these documents is incorporated herein by reference.

The above-described genome mapping strategies are based on mapping the distribution of short (from 6 bp to 8 bp) sequence motifs across the genome. However, the distribution of the sequence motifs is uneven at different genomic regions. Often, there are no appropriate sequence motifs in repetitive genomic regions, which results in large segments of the genome that cannot be mapped (Feuk, L., et al. (2006). “Structural variation in the human genome.” Nature Reviews Genetics 7(2): 85-97). Another challenge resides in detecting and typing structural variations or clinical diagnostics of specific structural variants. Target sequence-specific labeling of the structural variations is required to obtain accurate breaking points, but this cannot be achieved by sequence-motif mapping.

SUMMARY OF THE INVENTION

Provided herein are methods for sequence-specific labeling which is capable of targeting repetitive regions that often lack restriction site motif. In one aspect, methods of detecting a target nucleic acid sequence are provided. In one embodiment, the method includes contacting genomic DNA with a guide RNA having a portion complementary to the target sequence in the genomic DNA and with Cas9 nicking endonuclease (nickase) to produce a single-strand break (nick) in the genomic DNA at a specific location adjacent to the target sequence. The method further includes contacting the nicked DNA with a polymerase and fluorescently labeled nucleotide. The fluorescently labeled nucleotide is incorporated into the nicked DNA at the specific location and the target nucleic acid sequence is detected via fluorescent label.

In one embodiment, the method further includes contacting the genomic DNA with a second nicking endonuclease which is specific for a sequence motif in the genomic DNA thereby producing a second nick in the genomic DNA at the motif sequence. The nicked DNA is contacted with a polymerase and second fluorescently labeled nucleotide, wherein the second fluorescently labeled nucleotide is incorporated into the nicked DNA at the motif sequence location.

In another embodiment, the genomic DNA is contacted with a second nicking endonuclease which is specific for a sequence motif in the genomic DNA thereby producing a second nick in the genomic DNA at the motif sequence, wherein the first fluorescently labeled nucleotide is also incorporated into the nicked DNA at the motif specific location.

In another embodiment, the genomic DNA is contacted with multiple guide RNAs, each guide RNA having a portion complementary to a different target sequence in the genomic DNA, wherein each target nucleic acid sequence is detected via the same or different fluorescent label, thus providing a barcode of a portion of the genomic DNA.

Other aspects of the invention will be apparent from the description and examples below.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1A is a schematic graph illustrating sequence motif fluorescent DNA labeling. The repetitive elements contain no recognition motifs for nicking enzymes (empty circles), and therefore this region is undetectable by conventional sequence-motif dependent nick labeling. However, this region can be mapped by the new Cas9n fluorescent nick-labeling method (circles filled with lines).

FIG. 1B is a schematic graph illustrating target sequence specific fluorescent DNA labeling. The Cas9n fluorescent nick-labeling system can identify translocations by designing gRNAs to target Cas9n D10A to nick specific sequences near the breakpoints (diamonds). These translocation breakpoint regions can then be precisely mapped by visualizing the targeted Cas9n fluorescent nick-labeling system-dependent signals relative to the conventional sequence motif-label-dependent bar codes generated on the breakpoint-adjacent large DNA segments (circles).

FIG. 2 is a schematic graph illustrating Cas9n/gRNA target sequence specific fluorescent labeling for whole genome DNA mapping. The Cas9n fluorescent nick-labeling system uses a guide RNA (gRNA) to direct the Cas9 nuclease to a targeted site. The gRNA is composed of a trans-activating crRNA (tracrRNA) and a crRNA that contains a 20 nucleotide sequence that is complementary to the site of interest. A mutation in the RuvC-like domain nuclease alters the Cas9 enzyme to make only a single cut three nucleotides upstream of a protospacer adjacent motif (PAM) of the 3′-5′ strand of the target DNA. In the nick labeling method, after Cas9n D10A generates a nick, fluorophores are directly incorporated to the nick sites using Taq DNA Polymerase. These fluorophores can be detected using fluorescence microscopy.

FIG. 3A is a graph showing an example result of Cas9n-nick labeling. The expected labeling pattern and distances between labels is shown above corresponding fluorescent microscopy images. A single copy of HLS DUF1220 triplet is approximately 4.7 kb, which cannot be identified with the sequence-motif labeling method due to lack of recognition sites within the repeat region. The Cas9n fluorescent nick-labeling method was able to label this region using a gRNA designed to the DUF1220 triplet. These recognition sites were labeled with fluorescent labels. The tandem repeats were separated approximately 4.7 kb as predicted. The histogram of the DUF1220 gRNA labels is shown below the molecules. A total of 192 molecules were evaluated for the labeling efficiency.

FIG. 3B shows fluorescent microscopy images of three sgRNAs targeting three different locations (Gag, Pol, Env) within the HIV-1 genome. The sites were correctly labeled with the expected distances between each sgRNA. The image shows the linearized DNA after EcoRI digestion.

FIG. 3C shows fluorescent microscopy images of three gRNAs targeting the chromosome 1q subtelomere which generated the expected pattern.

FIG. 4A is a graph showing that the Cas9n fluorescent nick-labeling method was used to label the repetitive sequences of the DUF1220 triplet (circles filled with forward slashes). These molecules were then labeled with the conventional sequence motif labeling method with Nt.BspQI (circles filled with backslashes). The expected labeling pattern and distances between labels is shown above corresponding fluorescent microscopy images. The histogram of the labels is shown below the molecules. The red peaks represent the DUF gRNA labels and the green the Nt.BspQI labeling sites. A total of 43 molecules were evaluated for the labeling efficiency.

FIG. 4B shows a graph illustrating that the Cas9n fluorescent nick-labeling method was used to label the repetitive telomere sequences at the left end of these telomere-terminal DNA fragments. A gRNA targeting the repetitive telomere sequence (TTAGGG)_(n), correctly labeled the telomere region on the left end of this telomere-terminal 1q DNA fragment (circle filled with forward slashes). The sequence-motif labeling method with Nt.BspQI was performed after Cas9n/gRNA system labeling of telomeres, and is displayed with circle filled with backslashes.

FIG. 5A shows graphs illustrating Cas9n fluorescent nick-labeling created locus-specific and variant-specific barcodes. The expected labeling pattern and distances between labels are shown above corresponding fluorescent microscopy images. For each of the tested subtelomere regions specific gRNAs were designed and pooled in the in vitro nicking reaction mix in order to produce a unique pattern. In each instance, the expected pattern of the red fluorophores in the Cas9n fluorescent nick-labeling system was obtained.

FIG. 5B is a graph of human genomic DNA was labeled with Alu gRNAs, and followed by conventional sequence-motif labeling with Nt.BspQI. Conventional sequence-motif dependent labeling sites (black circles/ovals) were found infrequently within these genomic DNA regions. However, Alu-specific gRNAs directed Cas9n D10A to generate information-rich long-range barcoding patterns (white circles/ovals).

DETAILED DESCRIPTION OF THE INVENTION

The methods described herein used a modified CRISPR/Cas system. Clustered regularly interspaced short palindromic repeats (CRISPR) are segments of prokaryotic DNA containing short repetitions of base sequences. Each repetition is followed by short segments of “spacer DNA” from previous exposure to a bacteriophage virus or plasmid. CRISPR associated proteins (Cas) use the CRISPR spacers to recognize and cut these exogenous genetic elements in a manner analogous to RNA interference in eukaryotic organisms.

Recently, a new genome editing tool based on a bacterial CRISPR-associated protein-9 nuclease (Cas9) from Streptococcus pyogenes has been developed for generating double strand DNA breaks in vivo. To achieve site-specific DNA recognition and cleavage, the protein Cas9 must form a complex with a guide RNA (gRNA) comprised of a crRNA and a trans-activating crRNA (tracrRNA), which is partially complementary to the crRNA. The HNH and RuvC-like nuclease domains of Cas9 cut both DNA strands, generating double-stranded breaks (DSBs) at sites defined by a 20-nucleotide seed sequence within an associated crRNA transcript. Mutations of both nuclease domains generate nuclease-deficient Cas9 (dCas9) that is still capable of binding to gRNA and moving to the target sequence, and has been used to visualize repetitive DNA sequences. Other mutant forms which lack just the RuvC-like nuclease domain activity, only nick the DNA strand complementary to its crRNA, are is characterized as Cas9 nickases (Cas9n). This type of mutant of Cas9 has been used with paired singled guide RNA (sgRNA) targeting opposite strands of the same locus to generate DSBs with great precision.

As described herein, the Cas9n-dependent nicking protocol is combined with motif-dependent nicking to fluorescently label specific sequences for whole genome mapping and other applications as described herein, through in vitro nick-labeling.

Such Cas9n fluorescent nick labeling based sequence-specific labeling methods can be used to target repetitive regions which often lack appropriate restriction site motifs. This method can also help to precisely map the breaking points of structural variations such as translocations, by designing guide RNAs (gRNAs) to recognize and direct labeling of sequences near these breakpoints prior to high throughput single molecule analysis (FIG. 1B). Novel combinations of conventional sequence motif mapping and target-specific mapping will unleash the full potential of nanochannel array-based genome mapping methods.

Nicking

As used herein a “nickase” is an enzyme (e.g., an endonuclease) that causes breaks in one strand of the nucleic acid sequence (“nick”). Double-strand DNA breaks (DSB) occur or arise when both strands of the DNA duplex are severed. A “nick”, also known as single-strand DNA break (SSB), can stimulate gene correction without the problems of DSB repair because the uncut DNA strand acts as a template to permit healing without alteration of genetic material.

A Cas9 nickase (also called Cas9n or, alternatively, Cas9) is used in the methods herein. In one embodiment, the Cas9 nickase is Cas9 D10A. Cas9 H840A is another nickase useful in the methods described herein. In one embodiment, the Cas9 nickase is a mutant protein, which contacts the dsDNA and makes a nick of a single strand. The Cas9 nickase generates a single-strand DNA break at a specific location based on a guide RNA-defined target sequence, rather than a double-strand DNA break (“cut”) produced by the wild type enzyme. As used herein, the term “Cas9” is sometimes used interchangeably with “Cas9 nickase” or “Cas9 D10A”. The use of Cas9 nickase for targeting and labeling DNA has been described. See, e.g., WO 2016/028843, which is incorporated herein by reference.

It is contemplated that other nucleases engineered to create nicks can be utilized in place of, or in conjunction with, the Cas9 described herein. Such endonucleases include, without limitation, homing endonucleases (HE), meganucleases, Transcription activator-like effector nuclease (TALEN), Zinc finger nuclease (ZFN), prokaryotic Argonaute (pAgo), and BurrH-based nuclease (BuDN). Homing endonucleases (HE) are double stranded DNases that have large, asymmetric recognition sites (12-40 base pairs) and coding sequences that are extremely rare and usually embedded in either introns or inteins. Single base changes do not abolish cleavage by HE but reduce its efficiency to variable extents. As a result, the observed sequence specificity of HE is typically in the range of 10-12 base pairs.

Meganuclease are endodeoxyribonucleases characterized by a large recognition site (double-stranded DNA sequences of 12 to 40 base pairs), for example, I-SceI. Transcription activator-like effector nucleases (TALEN) are restriction enzymes that can be engineered to cut specific sequences of DNA. They are made by fusing a TAL effector DNA-binding domain to a DNA cleavage domain (a nuclease which cuts DNA strands). When combined with a nuclease, DNA can be cut at specific location. The restriction enzymes can be introduced into cells, for use in gene editing or for genome editing in situ.

Zinc-finger nucleases (ZFNs) are artificial restriction enzymes generated by fusing a zinc finger DNA-binding domain to a DNA-cleavage domain. Zinc finger domains can be engineered to target specific desired DNA sequences and this enables zinc-finger nucleases to target unique sequences within complex genomes. By taking advantage of endogenous DNA repair machinery, these reagents can be used to precisely alter the genomes of higher organisms and serve as a prominent tool in the field of genome editing.

Prokaryotic Argonaute (pAgo) is an endo-ribonuclease that uses a small RNA guide molecule to specifically target a complementary RNA transcript. Prokaryotic Argonautes are prokaryotic homologs of eukaryotic Argonaute proteins, which are key enzymes in RNA interference pathways. An Argonaute can bind and cleave a target nucleic acid by forming a complex with a designed nucleic acid-targeting nucleic acid. Cleavage can introduce double stranded breaks in the target nucleic acid. A nucleic acid can be repaired e.g. by endogenous non-homologous end joining (NHEJ) machinery. A piece of nucleic acid can be inserted. Engineering of non-genomic nucleic acid is also contemplated. Modifications of designed nucleic acid-targeting nucleic acids and Argonautes can introduce new functions to be used for genome engineering.

The method also includes contacting the genomic DNA with a motif-specific nicking endonuclease (also called “motif-specific nickase” or “motif-specific endonuclease”) thereby producing a second nick in the genomic DNA at the motif sequence. As used herein a “motif” sequence refers to a short DNA sequence, which generally recurs in the genome. One of the major functions of motifs is indicating sequence-specific binding sites for proteins such as nucleases and transcription factors (TF).

The motif-specific nickase is used to nick the DNA at various locations where the motif is present. Such motif-specific endonucleases are known in the art and are characterized in that they only cut one strand (“nick”) of the double stranded DNA, and are thus termed “nickases”. In one embodiment, the motif-specific nickase is Nt.BspQI. In another embodiment, the nickase is selected from Nt. CviPII, Nt.AlwI, Nt.BsmAI, Nt.BstNBI, Nb.BsmI, Nb.BsrDI, Nb.BtsI, Nb.BbvCI and Nt.BbvCI. Other motif-specific nickases useful herein are contained in Table 1 below, with the motif sequence that they recognize. The appropriate nickase may be selected based on desired application and the portion of the DNA being surveyed.

TABLE 1 Nickases NICKASE RECOGNITION SEQUENCE Nt.AlwI GGATC (4/none) Nb.BbvCI CCTCAGC (none/−2) Nt.BbvCI CCTCAGC (−5/none) Nt.BhaIII GAGTC (4/none) Nb.Bpu10I CCTNAGC (none/−2) Nt.Bpu10I CCTNAGC (−5/none) Nt.Bpu10IB CCTNAGC (−5/none) Nb.BsaI GGTCTC Nt.BsaI GGTCTC (1/none) Nb.BsmI GAATGC (none/−1) Nb.BsmAI GTCTC (none/5) Nt.BsmAI GTCTC (1/none) Nt.BsmBI CGTCTC (1/none) Nt.BspD6I GAGTC (4/none) Nb.BspQI GCTCTTC Nt.BspQI GCTCTTC (1/none) Nb.BsrDI GCAATG (none/0) Nb.BssSI CACGAG (none/−1) Nt.Bst9I GAGTC (4/none) Nt.BstNBI GAGTC (4/none) Nt.BstSEI GAGTC (4/none) Nb.BtsI GCAGTG (none/0) Nb.BtsCI GGATG Nt.BtsCI GGATG Nt.CviPII CCD (−3/none) Nt.CviQII RAG (−2/none) Nt.CviQXI RAG (−2/none) V.Gel16401III CGCG Nt.MlyI GAGTC (5/none) Nb.Mva1269I GAATGC (none/−1) Nb.SapI GCTCTTC Nt.SapI GCTCTTC

In one embodiment of the method, the genomic DNA is contacted with a guide RNA having a portion complementary to a target sequence in the genomic DNA and with Cas9 nickase to produce a nick in the genomic DNA at a specific location adjacent to the target sequence. The nicked DNA is contacted with a polymerase and fluorescently labeled nucleotide, wherein the fluorescently labeled nucleotide is incorporated into the nicked DNA at the specific location. The method further includes contacting the genomic DNA with a motif-specific nicking endonuclease thereby producing a second nick in the genomic DNA at the motif sequence. The nicked DNA is contacted with a polymerase and second fluorescently labeled nucleotide of the same color or different color, wherein the second fluorescently labeled nucleotide is incorporated into the nicked DNA at the motif sequence location. In another embodiment, the DNA is contacted with the motif-specific nickase prior to contacting with the Cas9 nickase. In another embodiment, where the same fluorescently labeled nucleotide is used, the DNA is contacted with both nickases at about the same time.

As used herein, the term “guide RNA” may refer to the combination of a long, chemically synthesized trans-activating CRISPR RNA (tracrRNA) plus a chemically synthesized CRISPR RNA (crRNA), which is partially complementary to the gene target site of interest. The term “guide RNA” may also, in another embodiment, refer to an expressed single guide RNA (sgRNA) that consists of both the crRNA and tracrRNA as a single construct.

In one embodiment, a guide RNA has a portion which is complementary to the target sequence in the genomic DNA. The “seed sequence” is complementary to the target DNA sequence and is incorporated into the crRNA. The crRNA is used in conjunction with the universal tracrRNA. In one embodiment, the gRNA is created by pre-incubating the tracrRNA and crRNA prior to contacting with the Cas9 enzyme.

In one embodiment, the crRNA and/or tracrRNA is chemically synthesized. In another embodiment, the crRNA and/or tracrRNA is in vitro transcribed. In yet another embodiment, the crRNA and/or tracrRNA is vector encoded and recombinantly produced.

In one embodiment, the gRNA (either sgRNA or crRNA/tracrRNA) is pre-incubated with the Cas9 before contacting with the genomic DNA to form a complex. In another embodiment, the gRNA and Cas9 are contacted with the DNA at approximately the same time.

In another embodiment, the guide RNA is expressed as a single guide RNA (sgRNA). In one embodiment, the sgRNA is chemically synthesized. In another embodiment, the sgRNA is in vitro transcribed. In yet another embodiment, the sgRNA is vector encoded and recombinantly produced.

Fluorescent Labeling

In the methods described herein, the DNA nicks are repaired using nucleotides which are fluorescently labelled, in conjunction with a DNA polymerase. The fluorescently labelled nucleotide is separately selected for both the Cas9 nick labeling step and the motif-specific nickase labeling step. That is, the same, or different fluorescent labelled nucleotides can be used for the Cas9 nick labeling step and the motif-specific nickase labeling step. Fluorescent labeling is a process of incorporating a fluorescent tag to a molecule or in a system to visualize the fluorescent tag, also known as a label or probe. Fluorescent dyes are covalently bound to biomolecules such as nucleic acids or proteins so that they can be visualized by fluorescence imaging. Suitable fluorescently labeled nucleotides are known in the art and include, without limitation, Alexa Fluor® 555-aha-dCTP, Alexa Fluor® 555-aha-dUTP, Alexa Fluor® 647-aha-dCTP, Alexa Fluor® 647-aha-dUTP, ChromaTide® Alexa Fluor® 488-5-dUTP, ChromaTide® Alexa Fluor® 546-14-dUTP, ChromaTide® Alexa Fluor® 568-5-dUTP, ChromaTide® Alexa Fluor® 594-5-dUTP, ChromaTide® Fluorescein-12-dUTP, ChromaTide® Texas Red®-12-dUTP, Fluorescein-aha-dUTP, DY-776-dNTP, DY-751-dNTP, ATTO 740-dNTP, ATTO 700-dNTP, ATTO 680-dNTP, ATTO 665-dNTP, ATTO 655-dNTP, OYSTER-656-dNTP, Cy5-dNTP, ATTO 647N-dNTP, ATTO 633-dNTP, ATTO Rho14-dNTP, ATTO 620-dNTP, DY-480XL-dNTP, ATTO 594-dNTP, ATTO Rho13-dNTP, ATTO 590-dNTP, ATTO Rho101-dNTP, Texas Red-dNTP, ATTO Thio12-dNTP, ATTO Rho12-dNTP, 6-ROX-dNTP, ATTO Rho11-dNTP, ATTO 565-dNTP, ATTO 550-dNTP, 5/6-TAMRA-dNTP, Cy3-dNTP, ATTO Rho6G-dNTP, DY-485XL-dNTP, ATTO 532-dNTP, 6-JOE-dNTP, ATTO 495-dNTP, BDP-FL-dNTP, ATTO 488-dNTP, 6-FAM-dNTP, 5-FAM-dNTP, ATTO 465-dNTP, ATTO 425-dNTP, ATTO 390-dNTP and MANT-dNTP. Suitable fluorescently labeled nucleotides also include dideoxynucleotides (ddNTPs). Each of the listed labels used with dNTPs is suitable for use with ddNTPs (e.g., ATTO 488-ddNTP) and is intended to refer to either a dNTP or ddNTP. Methods for nick-labeling are known in the art, and are described herein. See, e.g., Rigby, P. W. J., et al. [1977] J. Mol. Biol. 113:237, which is incorporated herein by reference.

In one embodiment, the fluorescent label used in the Cas9 labeling and the motif-specific labeling are different. In this embodiment, removal of the free nucleotides is performed. This can be accomplished by dialysis or via enzyme. Such methods are known in the art and are described in the examples. In one embodiment, the free nucleotides are removed via dialysis. In another embodiment, the free nucleotides are removed via enzyme. Such enzymes include shrimp alkaline phosphatase and pyrrolidonyl peptidase.

The nicked DNA is labeled using a DNA polymerase. In one embodiment, the DNA polymerase has exonuclease activity. Such DNA polymerases include, without limitation, Taq DNA Polymerase, E. coli DNA Polymerase I and Bst DNA Polymerase. In one embodiment, the polymerase is Taq DNA polymerase.

In one embodiment, the labeled DNA is repaired with a DNA ligase.

In another embodiment, the method includes RNAse treatment after the Cas9-dependent nicking step. In yet another embodiment, the method includes protease treatment.

In one embodiment, the DNA backbone is stained. Suitable dyes are known in the art and include, without limitation, YOYO-1, YOYO-3, YO-PRO-1, TOTO-1, TO-PRO-1, TO-PRO-3, TO-PRO-5, POPO, BOBO, JOJO, LOLO, ethidium bromide (EB), propidium iodide (PI), Hoechst 33342, 4′, 6-diamidino-2-phenylindole (DAPI), acridine orange, 7-AAD, LDS 751, hydroxystilbamidine, PicoGreen, OliGreen, RiboGreen, SYTOX Green/Blue/Orange, SYTO, and SYBR.

Visualizing DNA

In one embodiment, the genomic DNA is linearized and visualized after labeling. Methods of single DNA molecule fluorescence imaging are known in the art and include use of, e.g., nanochannels, nanopores or nanogaps.

In one embodiment, after the DNA is labelled using the methods described herein, the DNA is visualized. In one embodiment, the DNA labeling is visualized using a nanochannel, nanopore or nanogap. Visualizing DNA or optical mapping typically relies on sequence-specific DNA modifications at short target sites followed by imaging via Total Internal Reflection Fluorescence (TIRF) Microscopy (Chan, T. F., et al. (2006). “A simple DNA stretching method for fluorescence imaging of single DNA molecules.” Nucleic Acids Res 34(17): e1113), super-resolution microscopy and other microscopies. These techniques can roughly be divided into three groups: stretching over a surface, stretching via confinement in nanochannels, or stretching via elongational flow in micro/nanochannels. See, e.g., Noble, C., et al. (2015). “A fast and scalable kymograph alignment algorithm for nanochannel-based optical DNA mappings.” Plos One 10(4): e0121905); Chan, T. F., et al. (2006). “A simple DNA stretching method for fluorescence imaging of single DNA molecules.” Nucleic Acids Res 34(17): e113; Jing, J., et al. (1998). “Automated high resolution optical mapping using arrayed, fluid-fixed DNA molecules.” Proc Natl Acad Sci USA 95(14): 8046-8051, Das, S. K., et al. (2010). “Single molecule linear analysis of DNA in nanochannel labeled with sequence specific fluorescent probes.” Nucleic Acids Res 38(18): e177.; and Marie, R., et al. (2013). “Integrated view of genome structure and sequence of a single DNA molecule in a nanofluidic device.” Proc Natl Acad Sci USA 110(13): 4893-4898, each of which is incorporated herein by reference in its entirety.

Variants of flow stretching are implemented by attaching a large bead at the end of the DNA (Sriram, K. K., et al. (2014). “Direct optical mapping of transcription factor binding sites on field-stretched lambda-DNA in nanofluidic devices.” Nucleic Acids Res 42(10): e85) or by a DNA molecular tethering in both or a single end by chemical bonds (Qi, Z., et al. (2015). “DNA sequence alignment by microhomology sampling during homologous recombination.” Cell 160(5): 856-869), each of which is incorporated herein by reference in its entirety.

Nanopore and nanogap are also utilized for stretching, linearizing and imaging single DNA molecular. When a nanopore is immersed in a conducting fluid and a potential (voltage) is applied across it, an electric current due to conduction of ions through the nanopore can be observed. The amount of current is very sensitive to the size and shape of the nanopore. The capacitance, conductance and permittivity profiles of the Sub-10 nm nanogap electrodes are able to differentiate complementary, non-complementary and single mismatch target hybridization.

Elongation due to confinement in nanochannels can be performed by confining the DNA in one dimension by using a nanoslit or in two dimensions by using a nanoscale channel. See, e.g., Tegenfeldt, J. O., et al. (2004). “From the Cover: The dynamics of genomic-length DNA molecules in 100-nm channels.” Proc Natl Acad Sci USA 101(30): 10979-10983; and Lam, E. T., et al. (2012). “Genome mapping on nanochannel arrays for structural variation analysis and sequence assembly.” Nat Biotechnol 30(8): 771-776; Jo, K., et al. (2007). “A single-molecule barcoding system using nanoslits for DNA analysis.” Proc Natl Acad Sci USA 104(8): 2673-2678 and Kounovsky-Shafer, K. L., et al. (2013). “Presentation of large DNA molecules for analysis as nanoconfined dumbbells.” Macromolecules 46(20): 8356-8368, each of which is incorporated herein by reference.

Methods of using nanochannel arrays are known in the art. See, e.g., Lam et al, cited above). Briefly, a nanofluidic chip is provided that contains nanochannels that keep long DNA molecules in a consistent, uniformly elongated state. Fluorescently labeled DNA molecules are drawn into the nanochannels, held still and imaged automatically on the multicolor Irys® instrument. After imaging, additional sets of DNA molecules are streamed into the nanochannels for imaging. This process is repeated many times until the DNA is depleted or the nanochannels are rendered unusable as a result of clogging.

The nanofluidic chip contains three sets of nanochannels, each consisting of ˜4,000 channels that are 0.4 mm in length and 45 nm in diameter. Using 193-nm lithography in a nanofabrication process on the surface of a silicon substrate, nanochannel array chips are produced with precise diameters. DNA molecules in the 45-mm nanochannels cannot fold back on themselves and are forced by physical confinement to be in an elongated, linearized state.

As long DNA molecules in solution exist as coiled balls, a gradient region consisting of pillars and wider channels is placed in front of the nanochannels to allow the DNA molecules to uncoil as they flow toward the array. In this region, the physical confinement is sufficiently dense that the molecules are forced to interact with the pillars, yet sufficiently sparse that the DNA is free to uncoil. Once uncoiled, the DNA can then be efficiently flowed into the array in a linear manner, and visualized fluorescently.

As used herein, a “barcode” refers to a pattern of fluorescent labeling that is specific to a particular chromosome or portion of genomic DNA. For example, the specific pattern of fluorescence of the motif nick-labeling can identify the specific chromosome or portion of DNA being visualized. Alternatively, the pattern of nick labeling of both the motif-specific nick labeling and Cas9-dependent nick-labeling can be used to identity the specific chromosome or portion of DNA being visualized. Such barcodes are useful in the embodiment of the method in which the labeling is used to help assemble sequencing reads, e.g., for de novo sequencing, or to determine the location of integration of a viral sequence in the genome.

Telomere repeats are immediately distal to subtelomeric repeat elements (SREs) which are approximately 80% of the most distal 150 kb of human subtelomeres. Long SRE regions of about 150 kb have been identified in some alleles of various telomeres, whereas 7 telomeres have minimal or no SRE content. Read lengths of greater than 50 kb are required for assembling these regions using single-molecule sequencing, which is beyond the capability of current technology. Furthermore, the telomere repeat sequence (TTAGGG)n sequence lacks motif nicking sites recognized by currently available nicking endonucleases and therefore cannot be labeled with sequence-motif based methods.

As described herein, the CRISPR-Cas9 nick-labeling technology was combined with the nickase based nick-labeling procedure, to label telomeric repeats and subtelomeric regions separately. CRISPR-Cas9 labeling in conjunction with a global nickase enzyme motif-dependent labeling method allows for accurate telomere length measurement of each chromosome arm (FIG. 4B).

In the examples below, a mutant form of Cas9, termed Cas9 D10A, which is catalytically modified, was used with a telomere targeting gRNA to nick-label the telomeric repeats. The labeled DNA molecules were linearized in nanochannels and optically imaged. The labels from the enzyme motif-dependent labeling step were mapped to a reference to identify the chromosomal location of the corresponding Cas9 D10A labeled telomere. The intensity of the telomeric labeling was used and converted into base pairs to calculate the telomere length. This method was developed as a two-color and three-color approach, as described herein.

In one aspect, a method of detecting the length of an individual telomere is provided. As noted above, several embodiments are provided which utilize different schemes to produce different “barcodes” for identification the specific chromosome being identified. In one embodiment of the method, the first 20 nucleotides of the gRNA are complementary to the 3′ to 5′ telomere repeat region which is followed immediately by the protospacer adjacent motif (PAM) NGG.

In one embodiment, sometimes called the “3-color” method, the method includes contacting genomic DNA with a guide RNA having a portion complementary to the telomere repeat sequence in the genomic DNA and with a Cas9 nickase to produce a single-strand break (“nick”) in the genomic DNA at the telomere repeat sequence. The method further includes contacting the nicked DNA with a polymerase and fluorescently labeled nucleotide, wherein the fluorescently labeled nucleotide is incorporated into the nicked DNA at the telomere repeat sequence. The method further includes contacting the genomic DNA with a second nicking endonuclease which is specific for a sequence motif (“a motif-specific nickase”) in the genomic DNA thereby producing a second nick in the genomic DNA at the motif sequence. In one embodiment, the nick is made at the motif sequence in more than one location in the genome, e.g., throughout the genome. The nicked DNA is contacted with a polymerase and second fluorescently labeled nucleotide of different color, wherein the second fluorescently labeled nucleotide is incorporated into the nicked DNA at the motif sequence location. The length of the telomere is detected by measuring the fluorescence of first fluorescently labeled nucleotide at the telomere repeat location. The fluorescently labeled motif sequences are used as a barcode to identify the chromosome.

Methods

For each of the methods described herein, it is intended that the Cas9 nick-labeling may be performed before or after the motif-endonuclease nick labeling. Where multiple gRNAs are used in conjunction with Cas9 (e.g., for telomeric and subtelomeric sequences or multiple repeated sequences), these steps may be performed together or separate, with either step coming first.

In one embodiment, the method is useful in visualizing telomeres and determining the length thereof. Thus, in one embodiment, a guide RNA has a portion which is complementary to the telomere repeat sequence in the genomic DNA. There are 92 telomeres in a diploid human cell ranging between 0.5 kb and 20 kb, which are made up of the telomere repeat sequence (TTAGGG)n. In one embodiment, the seed sequence of 20 nucleotides complementary to the 3′-5′ strand of the telomere (UUAGGGUUAGGGUUAGGGUU—SEQ ID NO: 1) is incorporated into the crRNA. Other appropriate seed sequences based on the telomere repeat sequence can be designed by the person of skill in the art. The crRNA is used in conjunction with the universal tracrRNA. In one embodiment, the telomere gRNA is created by pre-incubating the tracrRNA and crRNA prior to contacting with the Cas9 enzyme. In one embodiment, the crRNA and/or tracrRNA is chemically synthesized. In another embodiment, the crRNA and/or tracrRNA is in vitro transcribed. In yet another embodiment, the crRNA and/or tracrRNA is vector encoded and recombinantly produced.

In one embodiment, the gRNA (either sgRNA or crRNA/tracrRNA) is pre-incubated with the Cas9 before contacting with the genomic DNA to form a complex. In another embodiment, the gRNA and Cas9 are contacted with the DNA at approximately the same time.

In another embodiment, the guide RNA is expressed as a single guide RNA (sgRNA). In one embodiment, the sgRNA is chemically synthesized. In another embodiment, the sgRNA is in vitro transcribed. In yet another embodiment, the sgRNA is vector encoded and recombinantly produced.

In another embodiment, sometimes called the “two-color” method, the genomic DNA is contacted with a guide RNA having a portion complementary to the telomere repeat sequence in the genomic DNA, and with a Cas9 nickase to produce a single-strand break in the genomic DNA at the telomere repeat sequence. The genomic DNA is also contacted with a second nicking endonuclease which is specific for a sequence motif in the genomic DNA, thereby producing a second nick in the genomic DNA at the motif sequence. The nicked DNA is contacted with a polymerase and fluorescently labeled nucleotide, wherein the fluorescently labeled nucleotide is incorporated into the nicked DNA at the telomere repeat sequence and at the motif sequence. The length of the telomere is detected by measuring the fluorescence of first fluorescently labeled nucleotide at the telomere repeat location and the fluorescently labeled motif sequences are used as a barcode to identify the chromosome. In this embodiment, the same fluorescent label is used to mark the motif sequences and the telomere sequences. In one embodiment, the telomere repeat sequence, the motif sequence, or both is present in more than one location in the genomic DNA.

In another embodiment, as an alternative, or in addition to, the motif-specific nickase, one or more additional guide RNAs are provided, which are complementary to one or more sequences in the genomic DNA other than the telomeric sequences. These gRNAs in combination with the Cas9 nickase can be used in a similar way to the motif-specific nickase to create a “barcode” to identify the chromosome. In one embodiment, the additional sequence is a sequence in the subtelomere region of the DNA. In another embodiment, multiple guide RNAs are provided to target multiple subtelomeric sequences. In one embodiment, the telomere repeat sequence, the subtelomeric sequence, or both is present in more than one location in the genomic DNA.

In one embodiment, the genomic DNA is contacted with a second guide RNA having a portion complementary to a sequence in the subtelomeric region of the genomic DNA and with Cas9 nickase to produce a single-strand break in the genomic DNA at the subtelomeric sequence. In one embodiment, the nicked DNA is contacted with a polymerase and second fluorescently labeled nucleotide, of a different color, wherein the second fluorescently labeled nucleotide is incorporated into the nicked DNA at the subtelomeric sequence location. The length of the telomere is detected by measuring the fluorescence of first fluorescently labeled nucleotide at the telomere repeat location, and the second fluorescently labeled subtelomeric sequence is used as a barcode to identify the chromosome.

In one embodiment, the second guide RNA comprises multiple guide RNAs, each guide RNA having a portion complementary to a different target sequence in the subtelomeric region of the DNA, wherein each subtelomeric sequence is detected via fluorescent label, thus providing a barcode of a portion of the genomic DNA.

In another embodiment, the DNA is nicked with both telomeric-directed Cas9 nickase and the subtelomeric-directed Cas9 nickase, and the same fluorescently labeled nucleotide is incorporated into the nicked DNA at the telomere repeat sequence and at the subtelomeric sequence. The length of the telomere is detected by measuring the fluorescence of the fluorescently labeled nucleotide at the telomere repeat location, and the fluorescently labeled subtelomeric sequence is used as a barcode to identify the chromosome.

As discussed herein, two labeling schemes for identifying genomic sequences are described, i.e., the so-called two-color and three-color methods. In both schemes, the same gRNA may be used to nick-label the first genomic sequence, and the remaining genomic DNA is globally nick-labeled using an endonuclease which is specific for a sequence motif in the DNA. Such motif-specific endonucleases are known in the art and are characterized in that they only cut one strand (“nick”) of the double stranded DNA, and are thus termed “nickases”. In one embodiment, the motif-specific nickase is Nt.BspQI. In another embodiment, the nickase is selected from Nt.CviPII, Nt.AlwI, Nt.BsmAI, Nt.BstNBI, Nb.BsmI, Nb.BsrDI, Nb.BtsI, Nb.BbvCI and Nt.BbvCI. Other motif-specific nickases useful herein are contained in Table 1 above.

In one embodiment, when two differently labeled fluorescent nucleotides are used, the unused nucleotides from the first nick labeling step are removed prior to the second nick labeling step.

In one embodiment, sometimes referred to as the “two-color” scheme, the Cas9n and Nt.BspQI are incubated with the human genomic DNA and simultaneously nick their respective regions. In the examples below, all of these nicks were labeled with green fluorescently labeled nucleotides using Taq DNA polymerase.

In the three color scheme, as described in the examples below, the Nt.BspQI was first incubated with the human genomic DNA and nicks the CGTCTTC motif. These nicks were labeled with green fluorescently labeled nucleotides using Taq DNA polymerase. After removing unused green fluorescent nucleotides with Shrimp Alkaline Phosphatase (SAP), the labeled DNA molecules were nicked with Cas9n and labeled with red fluorescently labeled nucleotides again using Taq DNA polymerase. The labeled DNA molecules were linearized in the nanochannels and optically imaged. De novo assembly was performed using the unique Nt.BspQI patterns referenced to hg38 and individual DNA molecules assembled are used for calculating the telomere length. This allowed for the identification of the chromosome arms with the corresponding telomere.

In one embodiment, the length of the telomere is detected by imaging the labeled telomere. In one embodiment, the contour of telomeric labeling is used to calculate telomere length. In another embodiment, the intensity of telomeric labeling is used to calculate telomere length. When viewed in a DNA imaging system, the longer the telomere, the more pixels it occupies. However, the ends of DNA molecules tend to fold back on to themselves which affects the length measurements. More importantly, due to photon scattering, even with single point emitter, several pixels collected photons. Thus, use of the contour method is less desirable for a telomere length less than 1 kb. This was observed with the telomeric labeling of the 8q fosmid with 800 bp of telomeric repeats occupying the same amount of pixels as a single fluorophore. Thus, for smaller telomeres, the total intensity of the telomere labels is a desirable method of measurement.

In one embodiment, to calculate the length of the telomere, the intensity of fluorescence is compared to a standard. For example, DNA which incorporates telomeres of known length can be included, visualized concurrently with the telomeres of unknown size, and the values compared. In another embodiment, the measurement is compared to a standard curve. Such values may be pre-existing based on prior data, or performed concurrently with the method. In another embodiment, background fluorescence is removed prior to calculation.

A decrease in the telomere length with later population doublings is detected because during each cell division telomeres become progressively shorter. It has been shown that a population of extremely short telomeres exists in human telomerase-positive human cancer cells and transformed fibroblast cells lacking tumor suppressor pathways. In fibroblasts with functional tumor suppressor pathways with or without telomerase these extremely short telomeres are rare.

Thus, in one aspect, a method of identifying a cancer cell is provided. The method includes detecting the length of the telomeres in the cell, or a portion thereof, using a method as described herein. The method further includes comparing the length of the telomeres to the length of telomeres to a standard. Such standards include a cancer cell, a healthy cell, an aging cell and a stem cell. In one embodiment, when one or more telomere is about 100 bp or more shorter as compared to a healthy cell, the cell is determined to be a cancer cell. In one embodiment, the method of identifying a cancer includes measuring the length of chromosome 8q.

The term “cancer” as used herein means any disease, condition, trait, genotype or phenotype characterized by unregulated cell growth or replication as is known in the art. A “cancer cell” is cell that divides and reproduces abnormally with uncontrolled growth. This cell can break away from the site of its origin (e.g., a tumor) and travel to other parts of the body and set up another site (e.g., another tumor), in a process referred to as metastasis. A “tumor” is an abnormal mass of tissue that results from excessive cell division that is uncontrolled and progressive, and is also referred to as a neoplasm. Tumors can be either benign (not cancerous) or malignant. The methods described herein are useful for the treatment of cancer and tumor cells, i.e., both malignant and benign tumors, so long as the cells to be treated have mitochondrial localization of the chaperones as described herein. In various embodiments of the methods and compositions described herein, the cancer can include, without limitation, breast cancer, lung cancer, prostate cancer, colorectal cancer, brain cancer, esophageal cancer, stomach cancer, bladder cancer, pancreatic cancer, cervical cancer, head and neck cancer, ovarian cancer, melanoma, acute and chronic lymphocytic and myelocytic leukemia, myeloma, Hodgkin's and non-Hodgkin's lymphoma, and multidrug resistant cancer. In one embodiment, the cancer is a drug resistant cancer.

In one embodiment, the methods described herein are used to enable mapping of the genome, or a portion thereof. In this embodiment, the gRNA is directed to a region which is repeated throughout the genome, or a portion thereof. The Cas9 is then used to nick-label these sequences in combination with a fluorescently labeled nucleotide and polymerase. This Cas9 nick-labeling is, in one embodiment, used in conjunction with motif nick-labeling to enable enhanced nanochannel array-based genome mapping.

Most resequencing projects rely on mapping the sequencing data to the reference sequence to identify variants of interest. When whole-genome assembly is attempted, paired-end sequencing of long DNA fragments provides scaffolds for assembly. As cloning of large DNA fragments is difficult, small-insert libraries of varying sizes may be prepared for paired-end sequencing, thus limiting the resolution of haplotypes and increasing the complexity, time and cost of the sequencing project. In addition, complex genomic loci, contain highly repetitive sequences and are particularly difficult to assemble. The methods described herein aid in de novo sequence assembly.

In one embodiment, a method of detecting a target nucleic acid sequence is provided. The method includes contacting genomic DNA with a guide RNA having a portion complementary to the target sequence in the genomic DNA and with Cas9 nickase to produce a single-strand break in the genomic DNA at a specific location adjacent to the target sequence. The method further includes contacting the nicked DNA with a polymerase and fluorescently labeled nucleotide, wherein the fluorescently labeled nucleotide is incorporated into the nicked DNA at the specific location. The target nucleic acid sequence is detected via fluorescent label.

In one embodiment, the guide RNA is directed to a sequence which is found in more than one location of the genome (a “repeated seqeuence”). In one embodiment, the guide RNA is directed to a SINE sequence. Short interspersed elements (SINEs) are highly repetitive sequences that retrotranspose into eukaryotic DNA through intermediates transcribed by RNA polymerase III (pol III). In many species, SINEs are a ubiquitously dispersed feature of the whole genome, often constituting a significant (˜10%) mass fraction of total DNA. SINEs cause mutations both by their retrotransposition within genes and by unequal recombination, and are widely considered to be examples of ‘selfish’ or ‘parasitic’ DNA. In one embodiment, the gRNA is directed to the Alu element. Alu sequences, with about one million copies, are the most abundant retroelements in humans, and account for up 10% of the human genome. These SINE sequences, exclusively found in primate genomes, have been particularly active in the human lineage even after human-chimpanzee divergence, where they are theorized to have contributed to some of the human-specific characteristics such as brain size (Britten, R. J. (2010). “Transposable element insertions have strongly affected human evolution.” Proceedings of the National Academy of Sciences of the United States of America 107(46): 19945-19948). In one embodiment, the gRNA sequence is directed to target a subset of the Alu sites. Other SINE sequences are known in the art and can be targeted using the methods described herein, depending on the subject genome being investigated, copy number, etc. A list of SINE sequences which can be targeted using the methods of the invention include those found in Table 2 below. Such sequences can be accessed at sines.eimb.ru/, which is incorporated herein by reference.

In another embodiment, the gRNA is directed to the DUF1220 protein domain. Genome sequences encoding DUF1220 protein domains have undergone an exceptional human lineage specific (HLS) increase in copy number, which was implicated in human brain size, pathology and evolution (See, Dumas, L. J., et al. (2012). “DUF1220-Domain Copy Number Implicated in Human Brain-Size Pathology and Evolution.” American Journal of Human Genetics 91(3): 444-454).

In one embodiment, the method includes the use of gRNAs to more than one repeated sequence. For example, one gRNA may be directed to the DUF1220 element and a second gRNA may be directed to the Alu element.

In another embodiment, the use of the Cas9 nick labeling of repeated sequences is used in combination with motif-specific nick labeling (e.g., thru motif-specific endonucleases).

In another aspect, the methods described herein are used to target viral sequences which have been incorporated into a host genome. For example, HIV-1 integrates into the human genome at different locations (See, Schroder, A. R. W., et al. (2002). “HIV-1 integration in the human genome favors active genes and local hotspots.” Cell 110(4): 521-529), each of which may be identifiable directly through an appropriate integration site labeling and global genome mapping strategy. Identification of such sites in HIV-1 latent cells is essential for understanding the molecular and epigenetic mechanisms underlying HIV-1 latency (See, Cohn, L. B., et al. (2015). “HIV-1 integration landscape during latent and active infection.” Cell 160(3): 420-432).

Similarly, this sequence-specific mapping approach is used to identify lentivirus random integration sites in the host cells, which may disrupt both endogenous gene expression and vector gene expression patterns (See, Desfarges, S. and A. Ciuffi (2010). “Retroviral integration site selection.” Viruses 2(1): 111-130).

As discussed in the Examples below, multiple sgRNAs were designed and tested to target HIV-1 structural region (Gag, Pol, Env) to determine the most effective gRNA that labels the HIV-1 genome. In one embodiment, the gRNA targets a structural region of the viral genome. In another embodiment, multiple gRNAs are used which target different regions of the viral genome. Other viruses which may be targeted include, without limitation, HIV, lentivirus, HPV, HSV, and HBV.

In another embodiment, the methods described herein are used to identify translocation breakpoint regions by designing gRNAs to target Cas9 nickase to nick specific sequences near the breakpoints (FIG. 1B). These translocation breakpoint regions can then be precisely mapped by visualizing the targeted Cas9n fluorescent nick-labeling system-dependent signals relative to the conventional sequence motif-label-dependent bar codes generated on the breakpoint-adjacent large DNA segments (light grey circles). Methods of determining the position of the breakpoint are known in the art and include maximum likelihood estimation. See, e.g., Chen et al, Mapping translocation breakpoints by next-generation sequencing, Genome Res. 2008 July; 18(7): 1143-1149, which is incorporated herein by reference. In one embodiment, the gRNAs are directed to specific sequences near the breakpoint locus.

The methods described herein can be modified to target any desirable genomic sequence. In such embodiments, the first set of gRNA together with Cas9 is designed to be complementary to, and target, any genomic sequence to generate a nick. The second set of gRNA-Cas9 complexes or nicking endonuclease can target any other desirable different genomic sequences than first set of gRNAs.

As used herein, the terms “complementary”, “complementarity” and “complement” have the same meaning as commonly used. Complementarity is a property shared between two nucleic acid sequences, such that when they are aligned antiparallel to each other, the nucleotide bases at each position will be complementary.

Technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention

TABLE 2 SINE Table—FIG. 6 SINE Length Σ: 213 (100%) (w/o tail) Taxon Distribution Copy number Mon1 270 Mammalia monotremes   2 × 10⁶ (platypus) Alu 282 Mammalia primates   1.1 × 10⁶ (human) 1.5 × 10⁵ (galago) CYN-II 154 Mammalia flying lemur   1 × 10⁶ (Cynocephalus) (with CYN-I & CYN-III) CYN-III 203 Mammalia flying lemur   1 × 10⁶ (Cynocephalus) (with CYN-I & CYN-II) PRE-1 246 Mammalia pigs and peccaries 10⁶ (pig) CYN-I 79 Mammalia flying lemur   1 × 10⁶ (t-SINE) (Cynocephalus) (with CYN-II & CYN-III) AfroSINE 230 Mammalia afrotherians   8 × 10⁵ (AfroSINEn) AfroLA 158 Mammalia elephants and mammoths   7 × 10⁵ (AFRO_LA, (Proboscidea) (mammoth) PSINE1) Lm1 200 Insecta Locusta, Schistocerca &   6 × 10⁵ (SGRP1) Chorthippus (locust) C 325 Mammalia Lagomorpha 5.4 × 10⁵ (CSINEn) Mar1 246 Mammalia marsupials   5 × 10⁵ (SINE1_MD) (Metatheria) (opossum) Cry 200 Reptilia Cryptodira   5 × 10⁵ (pol III/SINE, (straightneck turtles) Hirt) DANA 370 Actinopterygii zebrafish 4.5 × 10⁵ (Danio) Ther-1 252 Vertebrata mammals, birds & reptiles   4 × 10⁵ (MIR) (human)   1 × 10⁵ (mouse, rat) B4 278 Mammalia different rodent families   4 × 10⁵ (RSINE2, (mouse) ID_B1) 3.6 × 10⁵ (rat) Xt-1 208 Amphibia Xenopus tropicalis   4 × 10⁵ (MIR_Xt) ERE-B 235 Mammalia horse 3.8 × 10⁵ (ERE-2x, ERE- −176 (Equus caballus) 3x) B2 185 Mammalia Muridae, Cricetidae, 3.5 × 10⁵ (B2s, B2l, B3, Spalacidae & Rhizomiidae SINEB2) SINE type II 260 Mammalia bush babies and lorises   3 × 10⁵ (GarnAlu) (Lorisiformes) Lun1 291 Sarcopterygii lungfish   3 × 10⁵ (Lepidosiren paradoxa) Opo-1 193 Mammalia American marsupials 2.8 × 10⁵ (SINE-1_MD, SINE-2_MD) TUB 131 Mammalia tree shrews 2.5 × 10⁵ (TUB1, TUB2, (Scandentia) TUB3, TUB4) UnaSINE1 208 Actinopterygii eel 2.5 × 10⁵ (Anguilla japonica) Ped-2 185 Mammalia springhare 2.5 × 10⁵ (Pedetes capensis) DAS-III 387 Mammalia nine-banded armadillo 2.2 × 10⁵ (Dasypus novemcinctus) Tu type III 242 Mammalia tree shrews   2 × 10⁵ (TUBE1, TUS) (Scandentia) Bov-tA 204 Mammalia Bovidae   2 × 10⁵ (Bov-tAn) (cattle, goats & sheep) CAN 160 Mammalia Carnivora   2 × 10⁵ (SINEC*, MVB2) SINE type III 78 Mammalia bush babies and lorises   2 × 10⁵ (Garnel 1) (Lorisiformes) SS-1 121 Mammalia pig 1.9 × 10⁵ (SINE1_SS, (Sus scrofa) SUSINE2) ID-Spe 126 Mammalia Sciuromorpha 1.9 × 10⁵ (STRID2, STRID3, STRID4) Sepioth1 292 Cephalopoda Sepioteuthinae & 1.8 × 10⁵ Loliginidae Sepia 312 Cephalopoda cuttlefish 1.7 × 10⁵ (Sepiida) Bm1 430 Insecta silkworm 1.5 × 10⁵ (Bombyx mori) Sauria 350 Reptilia scaled reptiles 1.3 × 10⁵ (Squam1, ACA, (Squamata) AFE, POM, VIN) B1-dID 200 Mammalia Gliridae, Sciuridae & 1.1 × 10⁵ (STRIDE*) Aplodontidae (dormouse) 4.5 × 10⁴ (squirrel) 2.8 × 10⁴ (beaver) EM-1 126 Mammalia hedgehogs   1 × 10⁵ (Erine1) (Erinaceidae) TAL 237 Mammalia moles   1 × 10⁵ (SINE_UT) (Talpidae) vic-1 117 Mammalia Camelidae 10⁵ (ALPINE1, ALPINE2, SINE_VV) Talua 252 Insecta termites   1 × 10⁵ (Isoptera) DIP 190 Mammalia Dipodidae & Zapodidae   1 × 10⁵ (B2L_S) VES 191 Mammalia bat families 10⁵ (Ves2, Ves3, Vespertilionidae, Ves4) Molossidae, Phyllostomidae & Emballonuridae SOR 157 Mammalia shrews   1 × 10⁵ (Soricidae) MEN 259 Mammalia squirrels Menetes &   1 × 10⁵ Callosciurus ERI-2 186 Mammalia hedgehogs   1 × 10⁵ (Erine2) (Erinaceidae) WallSI4 233 Mammalia wallaby   1 × 10⁵ (Macropus eugenii) LF 453 Vertebrata coelacanths & tetrapods 10⁵ (LFSINE_Vert) (Sarcopterygii) (Latimeria menadoensis)   3 × 10² (mammals) IDL-Geo 190 Mammalia Geomyidae &   1 × 10⁵ (DIPODE1, Heteromyidae DIPODE2) Sepioth2B 229 Cephalopoda Sepioteuthinae & 9.4 × 10⁴ Loliginidae ERE-A 251 Mammalia horses   9 × 10⁴ (ERE-1x, ERE- (Equus spp.) 4x) Sepioth2A 293 Cephalopoda Sepioteuthinae   9 × 10⁴ OK 405 Cephalopoda Octopodidae 8.5 × 10⁴ Ther-2 220 Mammalia marsupials & placentals   8 × 10⁴ (MIR3) (human) Taluc 387 Insecta termites 7.9 × 10⁴ (Isoptera) DAS-IIb 355 Mammalia nine-banded armadillo 6.5 × 10⁴ (Dasypus novemcinctus) (with DAS-IIa) DAS-IIa 208 Mammalia nine-banded armadillo 6.5 × 10⁴ (Dasypus novemcinctus) (with DAS-IIb) CM-1 308 Chondrichthyes ghost shark 6.5 × 10⁴ (SINEX-1_CM, (Callorhinchus milii) SINE2-1_CM) Talud 261 Insecta termites 6.3 × 10⁴ (Isoptera) SP-D-Geo 225 Mammalia Geomyidae & 6.1 × 10⁴ (DIPODE3) Heteromyidae Ruka 230 Arachnida ixodid ticks   6 × 10⁴ (Rhipicephalus appendiculatus) Feilai 275 Insecta yellow fever mosquito   6 × 10⁴ (Aedes aegypti) Ped-1 231 Mammalia springhare   6 × 10⁴ (Pedetes capensis) Mar3 128 Mammalia marsupials 5.3 × 10⁴ (WSINE1) (Metatheria) (wallaby) TS 158 Angiospermae Solanaceae &   5 × 10⁴ Convolvulacea (tobacco) Mac1 311 Mammalia kangaroos 4.5 × 10⁴ (WALLSI2) (Macropodoidea) (wallaby) Cisc-1 290 Tunicata sea squirt   4 × 10⁴ (Ciona spp.) Fc-1 105 Mammalia Carnivora    4 × 10⁴- (tSINE_Fc)   6 × 10⁴ (cat, dog, panda) ACar-1 269 Reptilia green anole 3.9 × 10⁴ (SINE2- (Anolis carolinensis) 1_ACar) OR1 390 Cephalopoda Octopus spp. 3.3 × 10⁴ SlmI 227 Actinopterygii Salmonoidei   3 × 10⁴ RSINE1 139 Mammalia mouse & rat   3 × 10⁴ (mouse) 1.6 × 10⁴ (rat) DAS-I 76 Mammalia nine-banded armadillo 2.9 × 10⁴ (Dasypus novemcinctus) OR2 270 Cephalopoda Octopus spp. 2.8 × 10⁴ ID 75 Mammalia rodents  2.5 × 10⁴- (BC1, CAVID,   5 × 10⁵ MUSID, RDRE1_RN, RNALUIII, STRID1) Talub 277 Insecta termites   2 × 10⁴ (Isoptera) Bovc-Ta2 315 Mammalia Bovidae   2 × 10⁴ (cattle, goats & sheep) Anolis SINE2 456 Reptilia Green anole 1.4 × 10⁴ (Anolis carolinensis) MEG-RL 197 Mammalia megabats   1 × 10⁴ (Meg_PVa) (Megachiroptera) Au 180 Angiospermae Poaceae, Fabaceae,   1 × 10⁴ (SINE2-1_ZM) Solanaceae & Gramineae (Aegilops)   1 × 10² (other plants) CHRS 321 Mammalia cetaceans, hippopotamuses 10⁴-10⁵ (CHR-1, CHR- (120- ruminants & suiforms 2, CHRS-S, 340) SINE2-1_BT, CHRL, CHRL_BT, SINE2-1_BT, SINE2-2_BT, SINE2-3_BT) SlmII 465 Actinopterygii Salmonoidei   1 × 10⁴ SINE3 561 Amniota Amniota   1 × 10⁴ (AmnSINE1) (zebrafish) Oegop 371 Cephalopoda Oegopsida   1 × 10⁴ Mad-1 159 Angiospermae apple 9.5 × 10³ (SINE-1_Mad) (Malus × domestica) Pca-1 330 Mammalia hyrax 9.3 × 10³ (SINE2-1_Pca) (Procavia capensis) B1 135 Mammalia rodents    8 × 10³- (pB1, pB1dx, 6.5 × 10⁵ STR1) EC-1 398 Mammalia horse 7.7 × 10³ (SINE2-1_EC) (Equus caballus) Gecko 188 Insecta yellow fever mosquito 7.2 × 10³ (Aedes aegypti) MEG-RS 120 Mammalia megabats   7 × 10³ (Megachiroptera) FabaS-VII 124 Angiospermae barrel medic 6.9 × 10³ (Medicago truncatula) p-SINE1 115 Angiospermae rice   6 × 10³ (p-SINE2, p- (Oryza genus) (rice) SINE3, SINE1R5_OS, SINE1OS, SINE6_OS, SINE8_OS, SINE16_OS) STRIDM 199 Mammalia thirteen-lined ground 5.4 × 10³ squirrel (Spermophilus tridecemlineatus) BmSE 249 Insecta silkworm 5.2 × 10³ (Bombyx mori) Asu1 407 Insecta mosquito (Armigeres 5.1 × 10³ subalbatus) MEG-TR 182 Mammalia megabats   5 × 10³ (Megachiroptera) DR-2 174 Actinopterygii zebrafish   5 × 10³ (SINE2-1_DR, (Danio) CypS-I) DR-1 244 Actinopterygii zebrafish 4.8 × 10³ (SINE_DR1) (Danio) PTr-1 187 Angiospermae black cottonwood 4.1 × 10³ (SINE2-1_PTr, (Populus trichocarpa) SaliS-I) PTr-2 161 Angiospermae black cottonwood   4 × 10³ (SINE2-2_PTr, (Populus trichocarpa) SaliS-II) SP-7 199 Echinoidea Strongylocentrotus   4 × 10³ (SINE2-7_SP) purpuratus MEG-T2 212 Mammalia megabats   4 × 10³ (Megachiroptera) Lj-1 135 Angiospermae Lotus 3.1 × 10³ (Lj_SINE-1) japonicus & Medicago (Lotus japonicus) truncatula 3.9 × 10¹ (Medicago truncatula) CQ-6 208 Insecta mosquito 3.1 × 10³ (SINE6_CQ) (Culex pipiens) FabaS-IX 120 Angiospermae certain Fabales   3 × 10³ (Glycine max) 2.8 × 10² (Cajanus cajan &Medicago truncatula) FabaS-I 131 Angiospermae barrel medic 2.8 × 10³ (MT_SINE-1) (Medicago truncatula) EuphS-I 83 Angiospermae certain Euphorbiaceae 2.5 × 10³ (Manihot esculenta) 1.3 × 10² (Ricinus communis) CQ-3 404 Insecta mosquito 2.5 × 10³ (SINE3_CQ) (Culex pipiens) Lam1 230 Hyperoartia lamprey 2.4 × 10³ (Lethenteron reissneri) AFC-2 275 Actinopterygii African cichlids 2.4 × 10³ (SINE2-1_AFC) FabaS-II 116 Angiospermae barrel medic 2.2 × 10³ (SINE1_MT, (Medicago truncatula) MT-1) SP-4 262 Echinoidea Strongylocentrotus 2.1 × 10³ (SURF1, purpuratus SINE2-4_SP) Sine200 190 Insecta mosquito Anopheles   2 × 10³ (SINEX-1/2) gambiae SolS-I 179 Angiospermae Solanaceae   2 × 10³ (potato) 1.4 × 10³ (N. benthamiana) 6.3 × 10² (tomato) AFC 310 Actinopterygii African cichlids    2 × 10³- (SINE_AFC)   2 × 10⁴ SolS-II 216 Angiospermae Solanaceae 1.7 × 10³ (SINE2-1_STu) (potato, tomato) 2.2 × 10³ (N. benthamiana) FabaS-VIII 115 Angiospermae barrel medic 1.6 × 10³ (Medicago truncatula) HaSE2 284 Insecta certain owlet moths 1.6 × 10³ (HaSE2_DP) (Helicoverpa, Heliothis & (Danaus plexippus) Spodoptera) & monarch butterfly (Danaus plexippus) FabaS-VI 185 Angiospermae certain Fabales 1.4 × 10³ (Cajanus cajan) 4.6 × 10² (Medicago truncatula) CQ-2 238 Insecta mosquito 1.4 × 10³ (SINE2-2_CQ) (Culex pipiens) VitaS-I 93 Angiospermae grapevine 1.3 × 10³ (Vitis vinifera) SolS-IIIa 231 Angiospermae Solanaceae 1.3 × 10³ (potato)   1 × 10³ (tomato) SP-5 225 Echinoidea Strongylocentrotus 1.3 × 10³ (SPEC1, purpuratus SINE2-5_SP) SolS-V 107 Angiospermae Solanaceae 1.2 × 10³ (potato) 5.4 × 10² (tomato) 4.3 × 10² (N. benthamiana) OsSN2 283 Angiospermae rice & sorghum 1.2 × 10³ (SINE9_OS, (Oryza sativa & Sorghum (rice) SINE2-11_SBi, bicolor) SINE2-12_SBi) 2.2 × 10² (sorghum) Mad-2 193 Angiospermae apple   1 × 10³ (SINE-2_Mad, (Malus × domestica) SINE-4_Mad) ScroS-I 210 Angiospermae Common monkey-flower   1 × 10³ (Mimulus guttatus) Smα 332 Trematoda Schistosomatidae   1 × 10³ (Sjα, T2) HE1 340 Chondrichthyes higher elasmobranches 10³-10⁶ SolS-IV 193 Angiospermae Solanaceae   1 × 10³ (potato) 4.7 × 10² (N. benthamiana) 3.7 × 10² (tomato) SP-1 246 Echinoidea Strongylocentrotus 9.9 × 10² (SINE2-1_SP) purpuratus SolS-VI 225 Angiospermae Solanaceae 9.5 × 10² (N. benthamiana)   1 × 10² (potato) TguSINE1 124 Aves Estrildidae 9.4 × 10² CucuS-II 98 Angiospermae Cucurbitaceae 8.7 × 10² (Cucumis sativus) CQ-5 237 Insecta mosquito 8.4 × 10² (SINE5_CQ) (Culex pipiens) AFC-3 215 Actinopterygii African cichlids   8 × 10² (SINE3_AFC) OsSN3 176 Angiospermae rice 7.1 × 10² (CASINE) (Oryza sativa) CucuS-I 206 Angiospermae Cucurbitaceae 7.1 × 10² (Cucumis sativus) SaliS-V 268 Angiospermae black cottonwood   7 × 10² (Populus trichocarpa) FabaS-IV 205 Angiospermae barrel medic 6.9 × 10² (Medicago truncatula) FabaS-III 166 Angiospermae barrel medic 6.9 × 10² (MT_SINE-3) (Medicago truncatula) Xt-2 318 Amphibia Xenopus tropicalis   6 × 10² (SINE2-1_XT) OsSN1 285 Angiospermae rice 5.7 × 10² (F524) (Oryza sativa) Squam-2 191 Reptilia scaled reptiles 5.6 × 10² (Squamata) (anole) CQ-1 662 Insecta mosquito 5.5 × 10² (SINE2-1_CQ) (Culex pipiens) Lj-3 182 Angiospermae Lotus japonicus 5.5 × 10² (Lj_SINE-3) Twin 217 Insecta mosquito   5 × 10² (Culex spp.) MamSINE1 286 Mammalia mammals   5 × 10² (platypus) SaliS-III 168 Angiospermae black cottonwood 4.9 × 10² (Populus trichocarpa) BraS-I 134 Angiospermae Brassicaceae 4.6 × 10² (Arabidopsis lyrata) 3.2 × 10¹ (Arabidopsis thaliana) SB10 159 Angiospermae Brassica oleracea 455 (Bos d, e) SP-6 342 Echinoidea Strongylocentrotus 4.5 × 10² (SINE2-6_SP) purpuratus SB9 212 Angiospermae Brassica oleracea 395 (Bos c, f) GA-1 294 Actinopterygii stickleback 3.4 × 10² (SINE2-1_GA) (Gasterosteus aculeatus) Lj-2 191 Angiospermae Lotus japonicus 3.3 × 10² (Lj_SINE-2) SB1 171 Angiospermae Brassica oleracea 330 (S1) SaliS-IV 185 Angiospermae black cottonwood 3.2 × 10² (Populus trichocarpa) SolS-VII 211 Angiospermae Solanaceae   3 × 10² (N. benthamiana) 3.4 × 10¹ (tomato) FR-1 273 Actinopterygii fugu   3 × 10² (SINE_FR1) (Takifugu rubripes) SP-2 231 Echinoidea Strongylocentrotus 2.9 × 10² (SINE2-2_SP) purpuratus SB14 156 Angiospermae Brassica oleracea 285 (Bos j, k) PinS-I 227 Pinophyta spruces 2.7 × 10² (Picea spp.) PTr-3 183 Angiospermae black cottonwood 2.6 × 10² (Populus trichocarpa) FabaS-V 185 Angiospermae certain Fabales 2.4 × 10² (Medicago truncatula)  ~4 × 10¹ (Glycine max & Lotus japonicus) SB8 95 Angiospermae Brassica oleracea 215 Ras1 281 Actinopterygii Rasbora pauciperforata   2 × 10² CELE45 250 Nematoda Caenorhabditis elegans   2 × 10² AC1 380 Actinopterygii fugu & medaka   2 × 10² (SINE_FR2, (fugu) SINE_TE)   3 × 10³ (medaka) SolS-IIIb 246 Angiospermae Solanaceae 1.6 × 10² (potato, N. benthamiana)   1 × 10² (tomato) 5S-Sauria 236 Reptilia green anole 1.6 × 10² (Anolis carolinensis) SB12 170 Angiospermae Brassica oleracea 150 (Bos i1) SB2 149 Angiospermae Arabidopsis thaliana & 144 (RAthE1, SN2) Brassica oleracea (Arabidopsis) 165 (Brassica) PoaS-II 156 Angiospermae purple false brome 1.4 × 10² (Brachypodium distachyon) PoaS-I 156 Angiospermae purple false brome 1.4 × 10² (Brachypodium distachyon) CQ-4 583 Insecta mosquito 1.4 × 10² (SINE4_CQ) (Culex pipiens) SP-3 353 Echinoidea Strongylocentrotus 1.4 × 10² (SINE2-3_SP) purpuratus Sma I 150 Actinopterygii Salmonoidei    1 × 10²- (Fok I, pol   3 × 10⁴ III/SINE, SMAIDIV, SINE_SM) Cre-1 520 Chlorophyta Chlamydomonas   1 × 10² (SINEX-1_CR) reinhardtii Hpa I 200 Actinopterygii Salmonoidei 10²-10⁴ (OM2) SBi-1 282 Angiospermae Sorghum bicolor   1 × 10² (SINE2-1_SBi) MARE3 172 Mammalia mammals   1 × 10² (human) 1.4 × 10³ (opossum) Tu type II 278 Mammalia tree shrews 102 (TUBE4) (Scandentia) SB11 170 Angiospermae Brassica oleracea 100 (Bos G) Ava III 260 Actinopterygii Salmonoidei   1 × 10² UCON3 168 Euteleostomi Euteleostomi 102 Tu type I 179 Mammalia tree shrews 102 (TUBE2, (Scandentia) TUBE3, TUBE5) NV-1 257 Anthozoa starlet sea anemone 95 (SINE2-1_NV) (Nematostella vectensis) NV-2 324 Anthozoa starlet sea anemone 70 (SINE2-2_NV) (Nematostella vectensis) SB4 159 Angiospermae Arabidopsis thaliana 64 (RAthE3, SN1) NV-3 376 Anthozoa starlet sea anemone 60 (SINE2-3_NV) (Nematostella vectensis) SB3 297 Angiospermae Arabidopsis thaliana & 58 (RAthE2, Brassica oleracea (Arabidopsis) RBolE2) 510 (Brassica) SB13 225 Angiospermae Brassica oleracea 55 (Bos i2) SB15 206 Angiospermae Brassica oleracea 40 (Bos 1) SB6 297 Angiospermae Arabidopsis thaliana & 32 (ATSINE4) Brassica oleracea (Arabidopsis) 735 (Brassica) Cre-4 500 Chlorophyta Chlamydomonas 30 (SINEX-4_CR) reinhardtii Cre-3 496 Chlorophyta Chlamydomonas 30 (SINEX-3_CR) reinhardtii SB5 162 Angiospermae Arabidopsis thaliana & 7 (Bos a, ab, b, ai) Brassica oleracea (Arabidopsis) 645 (Brassica) Asu2 313 Insecta mosquito (Armigeres   6 × 10³ subalbatus) SB7 352 Angiospermae Arabidopsis thaliana & 3 Brassica oleracea (Arabidopsis) 210 (Brassica) EbuSINE1 345 Myxini hagfish (Eptatretus burgeri) Bf1SINE1 363 Leptocardii lancelet SacSINE1 454 Chondrichthyes dogfish shark (Squalus acanthias) Rhin-1 190 Mammalia bat families Rhinolophidae & Hipposideridae OS 469 Actinopterygii Salmonoidei SK 241 Cephalopoda squid (Loligo bleekeri) LmeSINE1 314 Sarcopterygii coelacanth (Latimeria menadoensis) Idio2 425 Cephalopoda Japanese pygmy squid (Idiosepius paradoxus) Idio1 259 Cephalopoda Japanese pygmy squid (Idiosepius paradoxus) EbuSINE2 345 Myxini hagfish (Eptatretus burgeri) MyrSINE 83 Mammalia anteaters (Myrmecophagidae) HaSE1 375 Insecta certain owlet moths (HzIS1) (Helicoverpa & Heliothis) Otri 252 Insecta mosquito Ochlerotatus triseriatus NymS-I 131 Angiospermae spatterdock (Nuphar advena) HaSE3 354 Insecta certain owlet moths (Helicoverpa, Heliothis, Trichoplusia & Mamestra) belongs and by reference to published texts, which provide one skilled in the art with a general guide to many of the terms used in the present application. The definitions contained in this specification are provided for clarity in describing the components and compositions herein and are not intended to limit the claimed invention.

It is to be noted that the term “a” or “an” refers to one or more. As such, the terms “a” (or “an”), “one or more,” and “at least one” are used interchangeably herein. The words “comprise”, “comprises”, and “comprising” are to be interpreted inclusively rather than exclusively. The words “consist”, “consisting”, and its variants, are to be interpreted exclusively, rather than inclusively. While various embodiments in the specification are presented using “comprising” language, under other circumstances, a related embodiment is also intended to be interpreted and described using “consisting of” or “consisting essentially of” language.

As used herein, the term “about” means a variability of 10% from the reference given, unless otherwise specified.

As used herein, genomic DNA refers to chromosomal DNA of the subject organism. Subject “organisms” or “subjects” or “patients” include, without limitation, an organism, such as an animal, such as Homo sapiens (human), Mus musculus (mouse), Rattus norvegicus (rat), Oryctolagus cuniculus (rabbit), Canis lupus familiaris (dog), Danio rerio (zebrafish), Drosophila melanogaster (fruit fly), Caenorhabditis elegans (roundworm), Sus scrofa (pig), Bos taurus (cow), Ovis aries (sheep), Capra aegagrus hircus (goat), non-human primates and insects; plant, such as crops and algae; fungi, bacterial, archaea, yeast, virus, protozoans, amebas and the like. Subject organisms can be selected by the person of skill in the art. It is also intended that genomic DNA can refer to a portion, rather than the totality, of the chromosomal DNA of an organism. Such portion may be contained within a plasmid, e.g. a bacterial artificial chromosome (BAC), or may be isolated.

The following examples illustrate several aspects and embodiments of the invention. As shown in the examples, the flexible and efficient fluorescent tagging of specific sequences allow us to obtain context specific sequence information along the long linear DNA molecules in the BioNano Genomics nanochannel. Our global nick-labeling scheme tags short recognition sequences, whose spatial relation can be translated into a genomic map. Not only can this integrated fluorescent DNA double strand labeling make the whole genome mapping more accurate, and provide more information, but it can also specifically target certain loci for clinical testing. Importantly, it renders the labeled double-stranded DNA available in long intact stretches for high-throughput analysis in nanochannel arrays as well as for lower throughput targeted analysis of labeled DNA regions using alternative methods for stretching and imaging the labeled large DNA molecules. Thus, this method will dramatically improve both automated high-throughput genome-wide mapping as well as targeted analyses of complex regions containing repetitive and structurally variant DNA.

We incorporated the CRISPR-Cas9 technology into our nick-labeling procedure for targeted, sequence-specific nicking and fluorescent labeling in one color, followed by global nickase enzyme motif-dependent labeling in a second color. We used two different forms of guide RNAs. For the HIV-plasmid experiments, the target sequence was cloned, amplified and then in vitro transcribed to result in an expressed single guide RNA (sgRNA). For all other experiments, the CRISPR RNAs (crRNAs) and transactivating CRISPRRNA (tracrRNA) were purchased from GE Dharmacon and pre-incubated to form each guide RNA (gRNA). In both cases, the gRNA-directed Cas9n D10A makes a precise single cut in one strand of the target double-strand DNA three nucleotides upstream of a protospacer adjacent motif (PAM) and fluorescent nucleotides are directly incorporated to these specific nick sites using Taq DNA Polymerase (FIG. 2). This approach promises to dramatically improve DNA mapping in complex and structurally variant genomic regions, as well as facilitate high-throughput automated whole-genome mapping.

We first established the Cas9n fluorescent nick-labeling conditions and investigated the labeling efficiency with BAC clones, fosmids and plasmids as model systems. FIG. 3A shows the Cas9n fluorescent nick-labeling results of HLS DUF1220 triplets on a BAC clone. The histogram of the label distribution is shown in the bottom graph of FIG. 3A. Genome sequences encoding DUF1220 protein domains have undergone an exceptional human lineage specific (HLS) increase in copy number, which is implicated in human brain size, pathology and evolution (22). A single copy of HLS DUF1220 triplet spans about 4.7 kb, and there are 12 copies on the BAC clone CH17-353B19. A gRNA probe was designed to target one unit of the triplets. Clearly, there are 12 copies detected on this 240 kb BAC clone and the distance between the each tripletmeasures 4.7 kb, which is in a good agreement with the clone sequence. The labeling efficiency of each locus was determined by evaluating 192 labeled BAC molecules, ranging from 87% to 98%. Interestingly, the middle copies (#6-8 labels from the left most label in FIG. 3A) are labeled at lower labeling efficiency. The Cas9n fluorescent nick-labeling is very specific.

The extra labels outside of the DUF domain were used to calculate the false positives in FIG. 3A. The spurious labels inside the DUF domain are harder to define because the DNA molecules inside the nanochannel are not static. Their slight movements and the limitation of optical resolution may cause inaccurate measurement of 4.7 kb repetitive sequences during the imaging time (200 ms). When all of the spurious spots inside and outside of the DUF domain were used, the false positive percentage is 0.6%. Next, we applied the Cas9n fluorescent nick-labeling method to a plasmid containing the HIV-1 genome. HIV-1 integrates into the human genome at different locations (23), each of which may be identifiable directly through an appropriate integration site labeling and global genome mapping strategy. Identification of such sites in HIV-1 latent cells is essential for understanding the molecular and epigenetic mechanisms underlying HIV-1 latency (24). Similarly, this sequence-specific mapping approach could be used to identify lentivirus random integration sites in the host cells, which may disrupt both endogenous gene expression and vector gene expression patterns (25). Multiple sgRNAs were designed and tested to target HIV-1 structural region (Gag, Pol, Env) to determine the most effective gRNA that labels the HIV-1 genome. The sites were correctly labeled with the expected distances between each sgRNA (FIG. 3B). However, the labeling efficiencies at each site of 36%, 58% and 44% from the left most label respectively, suggesting that labeling efficiency may be sequence- or region-dependent. It could also be the difference in labeling efficiency of sgRNA versus gRNA. The labeling efficiencies of the gRNAs were much higher than those achieved using sgRNAs. In a third model system, a fosmid containing a subtelomeric segment of human 1q ending in 100 bases of (TTAGGG)n was used to test the Cas9n fluorescent nick-labeling. We designed four guide RNA probes to target the (TTAGGG)n tract and three distinct loci on the subtelomere. The labeling pattern matches very closely with the positions of the gRNA seed sequences in the 1q reference sequence (FIG. 3C). However, the labeling efficiency is relatively lower for the telomere with 30%, while the labeling efficiency of subtelomeric markers are 95%, 79% and 99% respectively from the left most left label on 1q (FIG. 3C). This may be due to non-Watson-Crick pairings of the hexameric repeats, the high G+C content or the secondary structure of guide RNA probes targeting the (TTAGGG)n repeat.

We next tested the combination of Cas9n fluorescent nick-labeling with nicking endonuclease based sequence motif labeling. This approach has the potential to find wide applications in whole genome mapping of repetitive sequences as well as genotyping of structural variations and identification/mapping of viral integration sites. In FIG. 4A, the DUF1220 triplet repeats were first labeled with Cas9n fluorescent nick-labeling of red fluorescent nucleotides. These labeled DNA molecules were then globally nick-labeled with green nucleotides using Nt.BspQI to target the GCTCTTC motif. Twelve copies of DUF1220 triplets were detected spanning about 52 kb. Only the flanking regions of this 52 kb were shown to have the GCTCTTC motif, which can be mapped to reference genome to indicate the genomic locations of the DUF1220 triplet array. The histogram of the label distribution is shown in the bottom graph of FIG. 4A. Clearly, the combination of Cas9n fluorescent nick-labeling and nicking endonuclease sequence motif labeling, not only can detect the copy numbers of the DUF1220 triplets, but also can map the locations of the repeats.

The same approach was also applied to measure the telomere repeat length of a telomere-terminal fragment of chromosome 8q cloned in a fosmid. This fosmid carries 800 bp of the repetitive (TTAGGG)n sequence, which lacks motif nicking sites recognized by currently available nicking endonucleases and therefore cannot be labeled with current sequence-motif based methods. A gRNA specific for the telomere was designed. In FIG. 4B, the telomere was first labeled with Cas9n fluorescent nick-labeling of red fluorescent nucleotides. The labeled DNA molecules were then labeled with green nucleotides by Nt.BspQI nick-labeling to target the GCTCTTC motif. The telomere was correctly labeled on the end of the sequence. The length of this telomere region can be determined by measuring the length of the red fluorophore region and the intensity of its fluorescence relative to known controls after imaging. Sequence motif labeling over extended subtelomere regions linked on single large DNA molecules to the Cas9n fluorescent nick labeling (TTAGGG)n can be used to identify the specific subtelomere by comparison with genome-wide maps.

Cas9n fluorescent nick-labeling can be used to create locus-specific and variant-specific barcodes. We designed gRNAs to create barcodes to distinguish individual subtelomeres linked on single molecules to (TTAGGG)n tracts (FIG. 5A, Chr 1q and Chr 11q) and to distinguish variant copies of highly similar subtelomeric segmental duplications (Chr 15q and Chr 12p) using the Cas9n sequence specific labeling methods. It has been suggested that the shortest telomere or a small subset of the shortest telomeres in a cell determines the onset of senescence, apoptosis or genome instability (26, 27), the ability to systematically measure individual dysfunctionally short (TTAGGG)n tracts (rather than average telomere tract lengths in a sample as is typically measured currently) would provide important new high-resolution information on specific telomere functional status and identity at the single-molecule level. On the telomeric fosmids from Chr 15q and Chr 12p (FIG. 5A), Cas9n-gRNA directed nicks at two loci in a segmentally duplicated subtelomere region containing the WASH gene family (Linnardopoulou et al., 2007) were labeled in red. The same pair of gRNAs generated signals separated by 5.87 kb on the 15q fosmid and 7.5 kb on Chr 12p fosmid, as predicted by the reference sequences of these two highly similar but structurally non-identical regions. The patterns generated by three target-specific gRNAs on Chr1q and two gRNAs on 11 q are unique, easily distinguishable from each other and from both of the WASH gene-related patterns, and correspond exactly to what is expected from their respective reference sequences. These initial experiments indicate that specific gRNAs can be pooled to generate multiple specific nicks that, when labeled, can create custom barcodes predicted precisely by the respective reference sequences. FIG. 5B shows the results of experiments combining the Cas9n fluorescent nick-labeling and nicking endonuclease sequence motif labeling to map Alu elements in the human genome. Alu sequences, with about one million copies, are the most abundant retroelements in humans, and account for up 10% of the human genome. These SINE (Short Interspersed Nuclear Elements) sequences, exclusively found in primate genomes, have been particularly active in the human lineage even after human-chimpanzee divergence, where they are theorized to have contributed to some of the human-specific characteristics such as brain size.

A gRNA sequence was designed to target 280 000 Alu sites out of one million copies. Typical genomic DNA molecules are imaged and shown in FIG. 5B. One 180 kb molecule displays dense Alu elements with only two GCTCTTC motifs. Another DNA molecule shows dense Alu elements with two GCTCTTC motifs. The combined information can be used to map the DNA molecules to the reference genome and profile the distribution of Alu elements on the whole genome scale.

We have demonstrated the capabilities of our integrated approach of sequence motif fluorescent tagging and sequence specific labeling with Cas9n fluorescent nick-labeling system. The flexible and efficient fluorescent tagging of specific sequences allow us to obtain context specific sequence information along the long linear DNA molecules in the BioNano Genomics nanochannel. Our global nick-labeling scheme tags short recognition sequences, whose spatial relation can be translated into a genomic map. Not only can this integrated fluorescent DNA double strand labeling make the whole genome mapping more accurate, and provide more information, but it can also specifically target certain loci for clinical testing. Importantly, it renders the labeled double-stranded DNA available in long intact stretches for high-throughput analysis in nanochannel arrays as well as for lower throughput targeted analysis of labeled DNA regions using alternative methods for stretching and imaging the labeled large DNA molecules. Thus, this method will dramatically improve both automated high-throughput genome-wide mapping as well as targeted analyses of complex regions containing repetitive and structurally variant DNA. The publication McCaffrey et al, Nucleic Acids Research, 2015, 44(2):e11 doi: 10.1093/nar/gkv878 is specifically incorporated herein by reference in its entirety.

Example 1—Preparation of DNA Samples and Guide RNA

Target sequence-specific labeling with Cas9n fluorescent nick-labeling was carried out on the BAC clone CH17-353B19, fosmids carrying cloned telomere-terminal DNA fragments ending in several hundred bases of (TTAGGG)_(n)(Stong, N., et al. (2014). “Subtelomeric CTCF and cohesin binding site organization using improved subtelomere assemblies and a novel annotation pipeline.” Genome Res 24(6): 1039-1050), an HIV-1 entire genome-containing plasmid pEcoHIV-NL4-3-eLuc (gift from Dr Won-Bin Young at University of Pittsburgh) and genomic DNA isolated from human B-Lymphocyte cells NA12878 (Corriel Research Institute, NJ, USA).

The seed sequence of 20 nucleotides complementary to the 3′-5′ strand of the target template DNA were designed via a gRNA design tool (Feng Lab CRISPR Design Web Tool at crispr.mit.edu). Each seed sequence was incorporated into the crRNA. Two crRNAs for the genomic sequences of DUF1220 domain (in BAC clone), 1 for the telomere repeat sequence (TTAGGG)_(n) and 7 for subtelomeric sequences, along with the universal tracrRNA, were synthesized by GE Dharmacon. The fosmid and CH17-353B19 gRNAs were created by pre-incubating the tracrRNA (0.5 nmol) and corresponding crRNA (0.5 nmol) with 1×NEB Buffer 3 and 1×BSA at 4° C. for 30 min.

Three single guide RNAs (sgRNAs) containing seed sequence targeting HIV-1 structural gene regions (Gag, Pol and Env) were designed for efficiency and specificity using bioinformatics analysis tools. All the oligonucleotides for each target sequence (Table 3, Also see, Sequence Listing, which is incorporated by reference herein.) were synthesized in Alpha DNA (Montreal, Canada) and cloned into pKLVWG-sgRNA vector modified from pKLV-U6gRNA(BbsI)-PGKpuro2ABFP vector, a gift from Kosuke Yusa (Addgene plasmid #50946) (Koike-Yusa, H., et al. (2014). “Genome-wide recessive genetic screening in mammalian cells with a lentiviral CRISPR-guide RNA library.” Nat Biotechnol 32(3): 267-273). The vector was digested with BbsI and treated with Antarctic Phosphatase, and the linearized vector was purified with the QIAquick® nucleotide removal kit (QIAGEN®). A pair of oligonucleotides for each targeting site was annealed, phosphorylated and ligated to the linearized vector. The sgRNA expression cassette was validated by sequencing with U6 sequencing primer in GENEWIZ®. The validated vector was used as template for PCR with forward T7-U6 and reverse sgRNA primer to generate T7 promoter-driven gRNA expression cassette. Then the sgRNA for each target was in vitro transcribed using MEGAshortscript™ T7 transcription kit (Life Technology®). The quality of HIV-1 sgRNAs was verified by electrophoresis in 5% denaturing polyacrylamide gel.

TABLE 3 Sequences for gRNAs. Loci SEQ Loci number Sequence 5′-3′ Type ID NO. DUF1220 1-12 AAGUUCCUUUUAUGCAUUGG gRNA  2 HIV plasmid 1 CACCGCAGGATATGTAACTGACAG sgRNA  3 2 CACCGGCCAGATGAGAGAACCAAG sgRNA  4 3 CACCGAGAGTAAGTCTCTCAAGCGG sgRNA  5 Chr1q telomere 1 UUAGGGUUAGGGUUAGGGUU gRNA  6 Chr1q subtelomere 1 CCCCUGUUGCCAGAGCCAGU gRNA  7 2 GUAUUUAGUCAGAGGGCUAG gRNA  8 3 AUACAGUAGGAUAACCGCAA gRNA  9 Chr15q subtelomere 1 ACCUUGCUACCACGAGAGCA gRNA 10 2 UCCAUUGGUUUAAUUAGGAA gRNA 11 Chr11q subtelomere 1 GGUCCACCCUACAGAUGUGC gRNA 12 2 AGAUCAGCAGCCACGUGUGC gRNA 13 Chr12p subtelomere 1 ACCUUGCUACCACGAGAGCA gRNA 14 2 UCCAUUGGUUUAAUUAGGAA gRNA 15 Alu 1 UGUAAUCCCAGCACUUUGGG gRNA 16

Loci numbers designate the labels from left to right in FIGS. 3-5.

Example 2—Cas9n Fluorescent Nick-Labeling of Fosmids, HIV-1 Plasmid and BAC Clone CH17-353B19

The gRNAs or sgRNAs (5 μM) were incubated with 600 ng of Cas9n D10A (PNA Bio Inc), 1×NEB Buffer 3 and 1×BSA (NEB) at 37° C. for 15 min. The DNA (500 ng) was added to the mixture and incubated at 37° C. for 60 min. The nicked DNA was then labeled with 4.12 units of DNA Taq Polymerase (NEB), 0.1 μM of ATTO-532 dUTP dAGC and 1× Thermopol Buffer (NEB) at 72° C. 60 min. The labeled fosmids and BACs were cut and linearized with 5 units of NotI enzyme (NEB) at 37° C. for 60 min. The labeled pecoHIV-NL4-3-eLuc plasmid (17,099 bp) was digested with 20 units of a unique restriction enzyme EcoRI (at 5744 bp) (NEB). NotI and EcoRI were inactivated at 65° C. for 20 min.

The distances were calculated between spots using ImageJ. The histogram of the label distributions were plotted in Excel. If the pattern matched the predicted pattern we considered the labels as true positives. Missing labels were used for the calculation of labeling efficiency and the extra labels were used for calculating the false positive percentage.

To establish the Cas9n fluorescent nick-labeling conditions and investigated the labeling efficiency with BAC clones, fosmids and plasmids as model systems. The Cas9n fluorescent nick-labeling of HLS DUF1220 triplets on a BAC clone were assessed (FIG. 3A). A single copy of HLS DUF1220 triplet spanned about 4.7 kb, and there were 12 copies on the BAC clone CH17-353B19. A gRNA probe was designed to target one unit of the triplets. Clearly, there were 12 copies detected on this 240 kb BAC clone and the distance between the each triplet measured 4.7 kb, which was in a good agreement with the clone sequence. The labeling efficiency of each locus was determined by evaluating 192 labeled BAC molecules, ranging from 87% to 98%. Interestingly, the middle copies (#6-8 labels from the left most one in FIG. 3A) were labeled at lower labeling efficiency.

The Cas9n fluorescent nick-labeling was very specific. The extra labels outside of the DUF domain were used to calculate the false positives in FIG. 3A. The spurious labels inside the DUF domain were harder to define because the DNA molecules inside the nanochannel were not static. Their slight movements and the limitation of optical resolution may cause inaccurate measurement of 4.7 kb repetitive sequences during the imaging time (200 ms). When all of the spurious spots inside and outside of the DUF domain were used, the false positive percentage is 0.6%.

Further, the Cas9n fluorescent nick-labeling method was applied to a plasmid containing the HIV-1 genome. Multiple sgRNAs were designed and tested to target HIV-1 structural region (Gag, Pol, Env) to determine the most effective gRNA that labels the HIV-1 genome. The sites were correctly labeled with the expected distances between each sgRNA (FIG. 3B). However, the labeling efficiencies at each site of 36%, 58% and 44% from the left most label respectively, suggesting that labeling efficiency may be sequence- or region-dependent. It could also be the difference in labeling efficiency of sgRNA versus gRNA. The labeling efficiencies of the gRNAs were much higher than those achieved using sgRNAs.

In a third model system, a fosmid containing a subtelomeric segment of human 1q ending in 100 bases of (TTAGGG)_(n) was used to test the Cas9n fluorescent nick-labeling. Four guide RNA probes were designed to target the (TTAGGG)_(n) tract and three distinct loci on the subtelomere. The labeling pattern matched very closely with the positions of the gRNA seed sequences in the 1q reference sequence (FIG. 3C). However, the labeling efficiency was relatively lower for the telomere with 30%, while the labeling efficiency of subtelomeric markers were 95%, 79% and 99% respectively from the left most left label on 1q (FIG. 3C). This might be due to non-Watson-Crick pairings of the hexameric repeats, the high G+C content or the secondary structure of guide RNA probes targeting the (TTAGGG)_(n) repeat.

Example 3—The Two Color Genome Mapping with Cas9n Fluorescent Nick Labeling and Sequence-Motif Labeling

After nicking with Cas9n D10A as previously described in Example 2, the sample was digested with RNAseA (190 ng/μL, QIAGEN®) at 37° C. for 20 min. After digestion, the sample was labeled with ATTO 532-dATP, dTGC (100 nM) and 2.5 units of DNA Taq Polymerase (NEB) in the presence of 1× Thermopol Buffer (NEB) at 72° C. for 1 h. The sample was treated with 1 unit of SAP (USB® Products) and RNAseA (100 ng/μL) at 37° C. for 20 min and then 65° C. for 15 min. The nicks were repaired with 500 μM NAD+, 100 nM dNTPs and 20 kU of Taq DNA Ligase at 45° C. for 20 min. The sample was then treated with 6 mAU of QIAGEN Protease at 56° C. for 10 min and 70° C. for 15 min. The sample was dialyzed in TE on a 0.1 μm membrane (Millipore®) for 2 h. After dialysis, the sample was nicked with 10 units of Nt. BspQI (NEB) at 72° C. for 2 h. The nicked DNA was then labeled with 2.5 units of Taq DNA Polymerase (NEB), 0.1 μM ATTO-647 dUTP dAGC and 1× Thermopol Buffer (NEB) for 60 min at 72° C. The DNA backbone was stained with YOYO-1, and is shown in blue in all figures. The stained samples were loaded and imaged inside the nanochannels following the established protocol (See, Lam, E. T., et al. (2012). “Genome mapping on nanochannel arrays for structural variation analysis and sequence assembly.” Nature Biotechnology 30(8): 771-776). The distances were calculated between spots using ImageJ. The histogram of the label distributions were plotted in Excel. If the pattern matched the predicted pattern we considered the labels as true positives. Missing labels were used for the calculation of labeling efficiency and the extra labels were used for calculating the false positive percentage.

The combination of Cas9n fluorescent nick-labeling with nicking endonuclease based sequence motif labeling was tested. This approach was designed to find wide applications in whole genome mapping of repetitive sequences as well as genotyping of structural variations and identification/mapping of viral integration sites. In FIG. 4A, the DUF1220 triplet repeats were first labeled with Cas9n fluorescent nick-labeling of red fluorescent nucleotides. These labeled DNA molecules were then globally nick-labeled with green nucleotides using Nt.BspQI to target the GCTCTTC motif. Twelve copies of DUF1220 triplets were detected spanning about 52 kb. Only the flanking regions of this 52 kb were shown to have the GCTCTTC motif, which can be mapped to reference genome to indicate the genomic locations of the DUF1220 triplet array. The histogram of the label distribution was shown in the bottom graph of FIG. 4A. Clearly, the combination of Cas9n fluorescent nick-labeling and nicking endonuclease sequence motif labeling, not only can detect the copy numbers of the DUF1220 triplets, but also can map the locations of the repeats.

The same approach was also applied to measure the telomere repeat length of a telomere-terminal fragment of chromosome 8q cloned in a fosmid. This fosmid carries 800 bp of the repetitive (TTAGGG)_(n) sequence, which lacks motif nicking sites recognized by currently available nicking endonucleases and therefore cannot be labeled with current sequence-motif based methods. A gRNA specific for the telomere was designed. In FIG. 4B, the telomere was first labeled with Cas9n fluorescent nick-labeling of red fluorescent nucleotides. The labeled DNA molecules were then labeled with green nucleotides by Nt.BspQI nick-labeling to target the GCTCTTC motif. The telomere was correctly labeled on the end of the sequence. The length of this telomere region was determined by measuring the length of the red fluorophore region and the intensity of its fluorescence relative to known controls after imaging. Sequence motif labeling over extended subtelomere regions linked on single large DNA molecules to the Cas9n fluorescent nicklabeling (TTAGGG)_(n) can be used to identify the specific subtelomere by comparison with genome-wide maps.

Cas9n fluorescent nick-labeling was used to create locus-specific and variant-specific barcodes. gRNAs were designed to create barcodes to distinguish individual subtelomeres linked on single molecules to (TTAGGG)_(n) tracts (FIG. 5A, Chr 1q and Chr 11q) and to distinguish variant copies of highly similar subtelomeric segmental duplications (Chr 15q and Chr 12p) using the Cas9n sequence specific labeling methods. It has been suggested that the shortest telomere or a small subset of the shortest telomeres in a cell determines the onset of senescence, apoptosis or genome instability (Kaul, Z., et al. (2012). “Five dysfunctional telomeres predict onset of senescence in human cells.” Embo Reports 13(1): 52-59; and Zou, Y., et al. (2004). “Does a sentinel or a subset of short telomeres determine replicative senescence?” Molecular Biology of the Cell 15(8): 3709-3718), the ability to systematically measure individual dysfunctionally short (TTAGGG)_(n) tracts (rather than average telomere tract lengths in a sample as is typically measured currently) would provide important new high-resolution information on specific telomere functional status and identity at the single-molecule level. On the telomeric fosmids from Chr 15q and Chr 12p (FIG. 5A), Cas9n-gRNA directed nicks at two loci in a segmentally duplicated subtelomere region containing the WASH gene family (See, Linardopoulou, E. V., et al. (2007). “Human subtelomeric WASH genes encode a new subclass of the WASP family.” PLoS Genet 3(12): e237) were labeled in red. The same pair of gRNAs generated signals separated by 5.87 kb on the 15q fosmid and 7.5 kb on Chr 12p fosmid, as predicted by the reference sequences of these two highly similar but structurally non-identical regions. The patterns generated by three target-specific gRNAs on Chr1q and two gRNAs on 11q were unique, easily distinguishable from each other and from both of the WASH gene-related patterns, and corresponded exactly to what was expected from their respective reference sequences. These experiments indicate that specific gRNAs were pooled to generate multiple specific nicks that, when labeled, can create custom barcodes predicted precisely by the respective reference sequences.

FIG. 5B shows the results of experiments combining the Cas9n fluorescent nick-labeling and nicking endonuclease sequence motif labeling to map Alu elements in the human genome. A gRNA was designed to target 280 000 Alu sites out of one million copies. Typical genomic DNA molecules were imaged and shown in FIG. 5B. One 180 kb molecule displayed dense Alu elements with only two GCTCTTC motifs. Another DNA molecule showed dense Alu elements with two GCTCTTC motifs. The combined information can be used to map the DNA molecules to the reference genome and profile the distribution of Alu elements on the whole genome scale.

All publications cited in this specification are incorporated herein by reference in their entireties, as are U.S. Provisional Patent Application No. 62/410,322, filed Oct. 19, 2016, and U.S. Provisional Patent Application No. 62/410,324, filed Oct. 19, 2016. Similarly, the Sequence Listing filed herewith is hereby incorporated by reference. While the invention has been described with reference to particularly preferred embodiments, it will be appreciated that modifications can be made without departing from the spirit of the invention. Such modifications are intended to fall within the scope of the appended claims. 

1. A method of labeling a nucleic acid sequence comprising: a) contacting genomic DNA with a guide RNA having a portion complementary to a target sequence in the genomic DNA and with Cas9 nickase to produce a nick in the genomic DNA at a specific location adjacent to the target sequence; and contacting the nicked DNA with a polymerase and fluorescently labeled nucleotide, wherein the fluorescently labeled nucleotide is incorporated into the nicked DNA at the specific location; and b) contacting the genomic DNA with a motif-specific nicking endonuclease thereby producing a second nick in the genomic DNA at the motif sequence; and contacting the nicked DNA with a polymerase and second fluorescently labeled nucleotide of the same color or different color, wherein the second fluorescently labeled nucleotide is incorporated into the nicked DNA at the motif sequence location; wherein the nick-labeled DNA is detected via fluorescent label.
 2. The method of claim 1, wherein step b) is performed before step a).
 3. The method of claim 1, wherein the pattern of fluorescence provided by the Cas9 nickase is used as a barcode to identify the location or identity of the genomic DNA.
 4. The method of claim 1, wherein the pattern of fluorescence provided by the motif-specific nicking endonuclease is used as a barcode to identify the location or identity of the genomic DNA.
 5. The method of claim 1, wherein the pattern of fluorescence provided by the Cas9 nickase and the motif-specific nicking endonuclease is used as a barcode to identify the location or identity of the genomic DNA.
 6. The method according to claim 1, wherein the guide RNA comprises a crRNA and a tracrRNA.
 7. The method according to claim 1, wherein the guide RNA is a single gRNA sequence.
 8. The method according to claim 1, where the guide RNA and the Cas9 nickase are contacted with each other to form a complex, prior to contacting with the genomic DNA.
 9. The method according to claim 1, wherein the Cas9 nickase is Cas9 D10A or H840A.
 10. The method according to claim 1, wherein the genomic DNA is contacted with multiple guide RNAs, each guide RNA having a portion complementary to a different target sequence in the genomic DNA, wherein each target nucleic acid sequence is detected via the same or different fluorescent label, thus providing a barcode of a portion of the genomic DNA.
 11. The method of claim 1, wherein the guide RNA targets a repeated sequence.
 12. The method according to claim 1, wherein the target sequence is a viral nucleic acid sequence.
 13. The method according to claim 1, wherein the target sequence is a sequence adjacent to a translocation breakpoint.
 14. The method according to claim 1, further comprising ligating the labeled DNA with a ligase.
 15. The method according to claim 1, further comprising imaging the labeled genomic DNA in a nanochannel.
 16. The method according to claim 15, comprising transporting the genomic DNA into a nanochannel and maintaining the DNA in elongated form in the nanochannel.
 17. The method according to claim 1, further comprising staining the DNA backbone.
 18. The method according to claim 1, wherein the motif-specific nicking endonuclease is Nt.BspQI.
 19. The method according to claim 1, wherein the target sequence is present in more than one location in the genomic DNA.
 20. The method according to claim 1, further comprising one or more of: i. RNAse treatment; ii. Protease treatment; iii. Removing labeled free nucleotides; iv. Linearizing the DNA; and v. Fragmenting the DNA. 